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We study the processes e + e~ — > 2(7r + 7r~)7r°7, 2{-K + -K~)r)^, K + K~ '7r + ~tt~ V : '7 and K + K~-K + -K~ru 
with the hard photon radiated from the initial state. About 20000, 4300, 5500 and 375 fully 
reconstructed events, respectively, are selected from 232 fb _1 of BABAR data. The invariant mass 
of the hadronic final state defines the effective e + e~ center-of-mass energy, so that the obtained 
cross sections from the threshold to about 5 GeV can be compared with corresponding direct 
e + e~ measurements, currently available only for the rj-K + ir~ and lutt + tt~ submodes of the e + e _ — ► 
2(7T + 7r~)7r channel. Studying the structure of these events, we find contributions from a number 
of intermediate states, and we extract their cross sections where possible. In particular, we isolate 
the contribution from e + e~ — ► u;(782)7r + 7r~ and study the oj(1420) and w(1650) resonances. In the 
charmonium region, we observe the J/tp in all these final states and several intermediate states, as 
well as the ip(2S) in some modes, and we measure the corresponding branching fractions. 

PACS numbers: 13.66.Bc, 14.40.Cs, 13.25.Gv, 13.25.Jx, 13.20.Jf 
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I. INTRODUCTION 

Electron-positron annihilation at fixed center-of-mass 
(cm.) energies has long been a mainstay of research in 
elementary particle physics. The idea of utilizing initial- 
state radiation (ISR) to study e + e~ reactions below the 
nominal cm. energies was outlined in Ref. [1], and dis- 
cussed in the context of high-luminosity <p an d B fac- 
tories in Refs. [2-4]. At high energies, e + e~ annihi- 
lation is dominated by quark-level processes producing 
two or more hadronic jets. However, low-multiplicity 
processes dominate at energies below about 2 GeV, and 
the region near charm threshold, 3.0-4.5 GeV, features 
a number of resonances [5]. These allow us to probe a 
wealth of physics parameters, including cross sections, 
spectroscopy and form factors. 

Of particular current interest are several recently ob- 
served charmonium states and a possible discrepancy be- 
tween the measured value of the anomalous magnetic mo- 
ment of the muon (g — 2)^ and that predicted by the 
Standard Model [6] . Charmonium and other states with 
can be observed as resonances in the cross 
section, and intermediate states may be present in the 
hadronic system. Measurement of the decay modes and 
their branching fractions is important in understanding 
the nature of these states. The prediction for (g — 2)^ is 
based on hadronic-loop corrections measured from low- 
energy e + e~ — > hadrons data, and these dominate the 
uncertainty on the prediction. Improving this predic- 
tion requires not only more precise measurements, but 
also measurements over the entire energy range and in- 
clusion of all the important subprocesses in order to un- 
derstand possible acceptance effects. ISR events at B 
factories provide independent measurements of hadronic 
cross sections with complete coverage from the produc- 
tion threshold to about 5 GeV. 

The cross section for the radiation of a photon of 
energy E 1 followed by the production of a particular 
hadronic final state / is related to the corresponding di- 
rect e + e~ — > / cross section cr/(s) by 



da 1 f{s,x) 
dx 



= W{s,x)-a f {s(l-x)) 



(1) 



where s/s is the initial e + e~ cm. energy, x = 2E 1 j \fs 
is the frac tional energy of the ISR photon and £ c . m . = 
\/s(l — x) is the effective cm. energy at which the fi- 
nal state / is produced. The probability density function 
W(s,x) for ISR photon emission has been calculated with 
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better than 1% precision (see e.g. Ref. [4]). It falls rapidly 
as E-y increases from zero, but has a long tail, which com- 
bines with the increasing af(s(l— x)) to produce a sizable 
cross section at very low E c , m ,. The angular distribution 
of the ISR photon peaks along the beam directions, but 
10-15% [4] of the photons are within a typical detector 
acceptance. 

Experimentally, the measured invariant mass of the 
hadronic final state defines -E c m .. An important feature 
of ISR data is that a wide range of energies is scanned 
simultaneously in one experiment, so that no structure 
is missed and the relative normalization uncertainties in 
data from different experiments or accelerator parame- 
ters are avoided. Furthermore, for large values of x the 
hadronic system is collimated, reducing acceptance issues 
and allowing measurements at energies down to produc- 
tion threshold. The mass resolution is not as good as a 
typical beam energy spread used in direct measurements, 
but the resolution and absolute energy scale can be mon- 
itored by the width and mass of well known resonances, 
such as the J/tjj produced in the reaction e + e~ — > J/ipl- 
Backgrounds from e+e~ — > hadrons events at the nomi- 
nal y/s and from other ISR processes can be suppressed 
by a combination of particle identification and kinematic 
fitting techniques. Studies of e + e~ — > ^+^^7 and sev- 
eral multi-hadron ISR processes using BABAR data have 
been reported [7—11], demonstrating the viability of such 
measurements . 

The contributions to the 2(7r + 7r~)7r° final state from 
the r)Tr + TT~ and ujit + tt~ channels have been measured di- 
rectly by the DM1 [12], DM2 [13, 14], CMD2 [15] and 
ND [16] collaborations for \/s < 2.2 GeV. In this paper 
we present a comprehensive study of the 2(7r + 7r~)7r° fi- 
nal state along with new measurements of the 2(7r + 7r~)?7, 
K + K~ it + -k° and K + K^-K + -K^rj final states. In all 
cases we require detection of the ISR photon and per- 
form a set of kinematic fits. We are able to suppress 
backgrounds sufficiently to study these final states from 
their respective production thresholds up to 4.5 GeV. In 
addition to measuring the overall cross sections, we study 
the internal structure of the events and measure cross sec- 
tions for a number of intermediate states. We study the 
charmonium region, measuring several J/ip and ip(2S) 
branching fractions. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric energy e + e~ 
storage rings. The total integrated luminosity used is 
232 fb"\ which includes 211 ib -1 collected at the T(4S') 
peak, x/s = 10.58 GeV, and 21 ffcT 1 collected below the 
resonance, at \/s = 10.54 GeV. 

The BABAR detector is described elsewhere [17]. Here 
we use charged particles reconstructed in the tracking 
system, which compresed the five-layer silicon vertex 
tracker (SVT) and the 40-layer drift chamber (DCH) in 
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a 1.5 T axial magnetic field. Separation of charged pi- 
ons, kaons and protons uses a combination of Chcrcnkov 
angles measured in the detector of internally reflected 
Cherenkov light (DIRC) and specific ionization measured 
in the SVT and DCH. Here wc use a kaon identification 
algorithm that provides 90-95% efficiency, depending on 
momentum, and rejects pions and protons by factors of 
20-100. Photon and electron energies are measured in 
the CsI(Tl) electromagnetic calorimeter (EMC). 

To study the detector acceptance and efficiency, we use 
a simulation package developed for radiative processes. 
The simulation of signal and background hadronic fi- 
nal states is based on the approach suggested by Czyz 
and Kiihn [18]. Multiple soft-photon emission from 
the initial-state charged particles is implemented with 
a structure-function technique [19, 20], and photon ra- 
diation from the final-state particles is simulated by the 
PHOTOS package [21]. The accuracy of the radiative 
corrections is about 1%. 

We simulate the 2(7r + 7r~)7r final state both accord- 
ing to phase space and with models that include the 
r/p, LJTr + ir~ and ujf (980) intermediate states, and the 
K + K~ir + ir~ir° final state both according to phase space 
and including the intermediate <f> and/or r/ resonances. 
The generated events go through a detailed detector sim- 
ulation [22] , and we reconstruct them with the same soft- 
ware chain as the experimental data. Variations in de- 
tector and background conditions are taken into account. 

We generate a large number of background pro- 
cesses, including the ISR channels e + e~^7r + 7r~7r + 7r~7, 
K + K~7r + iT~-f, 2(ir + ir~)ir°ir°j, and K + K~ir + ir~ir w j. 
These can contribute due to a combination of particle 
misidentification, and missing or spurious tracks or pho- 
tons. In addition, we study the non-ISR backgrounds 
e + e~^qq (q = u, d, s, c) generated by JETSET [23] and 
e+er-fT+T- by KORALB [24]. The contribution from 
the T(4S) decays is negligible. The cross sections for 
these processes are known with about 10% accuracy or 
better, which is sufficient for these measurements. 



III. EVENT SELECTION AND KINEMATIC FIT 

In the initial selection of candidate events, we con- 
sider photon candidates in the EMC with energy above 
0.03 GeV and charged tracks reconstructed in the DCH 
or SVT or both that extrapolate within 0.25 cm of the 
beam axis and within 3 cm of the nominal collision point 
along the axis. We require a high-energy photon in the 
event with an energy in the initial e + e~ cm. frame of 
E 1 > 3 GeV, and exactly four charged tracks with zero 
net charge that combine with a pair of other photons 
to roughly balance the momentum of the highest-energy 
photon. We fit a vertex to the set of charged tracks and 
use it as the point of origin to calculate the photon direc- 
tions. Most events contain additional soft photons due 
to machine background or interactions in the detector 
material. 



We subject each of these candidate events to a set of 
constrained kinematic fits, and use the fit results, along 
with charged-particle identification, both to select the fi- 
nal states of interest and to measure backgrounds from 
other processes. We assume the photon with the high- 
est E 7 is the ISR photon, and the kinematic fits use its 
direction and energy along with the four-momenta and 
covariance matrices of the initial e + e~ and the set of 
selected tracks and photons. The fitted three- momenta 
for each charged track and photon arc used in further 
kinematical calculations. 

We pair all non-ISR photon candidates and consider 
combinations with invariant mass within ±30 McV/c 2 of 
the tt°(t]) mass as 7r°(?7) candidates. For each candidate 
event, we perform a kinematic fit using the momenta of 
the ISR photon, of the two additional photons and of the 
four tracks (under the relevant mass hypotheses for each 
considered final state), imposing five constraints (5C): 
the two photons invariant mass must match the nominal 
7T° or r\ mass, and the energy and momentum of the whole 
event must match the energy and momentum of the ini- 
tial e+e - state. We retain the combination with the low- 
est x 2 , either xl nn o or xl nv > as a 2(7r+7r~)7r° or 2("7r+7r~)?7 
candidate, respectively. If the four tracks include ex- 
actly one identified K + and one K~, we perform a set of 
similar fits under the K + K ~7r + 7r~7r° and K + K~-K + -K~r] 
hypotheses and retain the two-photon combination with 
the lowest X2K2inr /ri- To re duce the large background 
from ISR 2(tt + it~) events [9] with no additional neu- 
tral hadrons, we also fit each candidate event under the 
2(7r+7r~) hypothesis and require xl T >20. 

IV. THE 2(7T + 7r-)7r° FINAL STATE 

A. Final Selection and Backgrounds 

To suppress K + K ~ir + TT~ir° background, we require 
that no more than one track in the event is identi- 
fied as a kaon, and we also fit under all four possible 
K + K~ir + ir~iT hypotheses and require X2K2mr° > ^0- To 
reject any K ± K°ir T n° background, we require all tracks 
to extrapolate within 2.5 mm of the beam axis. The 
result of the 5C fit for the remaining events under the 
2(7T + 7r~)7r°7 hypothesis with the 2(7r + 7r~)7r invariant 
mass up to 4.5 GeV/c 2 is used for the final event selec- 
tion and background estimate. We consider two types of 
background: a non-ISR type background and ISR-type 
background. 

The non-ISR type background comes from the e + e~^ 
qq events and we estimate it using the JETSET simu- 
lation. It is dominated by events with a hard it pro- 
ducing a fake ISR photon, and the similar kinematics 
causes it to peak at low values of x\^o- We evaluate 
this background in a number of E c _ m , ranges by com- 
bining the ISR photon candidate with another photon 
candidate in both data and simulated events, and com- 
paring the 7T° signals in the resulting 77 invariant mass 
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20 40 60 80 100 

% 2 (47C7C°) 

FIG. 1: Distribution of \ from the five-constraint fit for 
2(7r + 7r - )7r° candidates in the data (points) after subtracting 
the e + e~ — > qq background (hatched histogram). The open 
histogram is the distribution for simulated signal events, nor- 
malized as described in the text. The dashed histogram is the 
estimated backgrounds from other ISR channels, as described 
in the text. 



distributions. The simulation gives an -E c .m. -dependence 
consistent with the data, so we normalize it by an overall 
factor. The hatched histogram in Fig. 1 represents this 
background and we subtract it from the experimental 
distribution. 

The xIttttO distribution for the remaining events is 
shown in Fig. 1 as points, and the open histogram is the 
distribution for the simulated 2(tt + tt^)tt° events. The 
simulated distribution is normalized to the data in the 
region xJLtt < 10 where the backgrounds and radiative 
corrections are smallest. The experimental distribution 
has contributions from ISR-type background processes, 
but the simulated distribution is also broader than the 
expected 5C x 2 distribution. This is due to multiple soft- 
photon emission from the initial state and radiation from 
the final-state charged particles, which are not taken into 
account by the fit. The shape of the \ 2 distribution at 
high values was studied in detail [9, 10] using ISR pro- 
cesses for which very clean samples can be obtained with- 
out any limit on the \ 2 value and MC signal events have 
been found to accurately simulate it. 

All ISR-type background sources are consistent with 
having a x\^a distribution that is nearly uniform over 
the range shown in Fig. 1. As an example, the Xa-~o 
distribution predicted from our simulations of other ISR 
channels (see Sec. II) is shown as the dashed histogram, 
with the main contribution from the 2(7r + 7r _ )7r°7r°7 pro- 
cess [10]. We therefore determine the xl^o distribution 
of the ISR-type background (x"jsRbkg) from the data dis- 
tribution, by subtracting the x\^^o distributions of simu- 



lated signal and of the qq backgrounds, both of which are 
normalized to data as mentioned above. The obtained 
Xjsmkg distribution is in agreement with simulation in 
shape, but contains events from the processes which are 
not included into simulation. 

In order to determine the mass spectrum of the gen- 
uine 2(7r + 7T~)7r° events, we define signal (x^- o <40) and 
control (40< xL-tt <80) regions as shown in Fig. 1. The 
signal region of Fig. 1 contains 30776 data and 17477 
simulated signal events, and the control region contains 
11829 data and 2012 simulated events. For each mass 
bin, the number of signal and ISR-background events in 
the signal region are extracted using the observed num- 
bers of events in the two regions with the qq background 
subtracted, and the two ratios of contributions expected 
from the shapes of the simulated signal and the XiSRbkg 
distributions. The qq subtraction is actually performed 
using a smooth function interpolating the simulated mass 
distribution. 
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m(2(7i + 7i")7i°) (GeV/c 2 ) 

FIG. 2: Invariant mass distribution for selected 2(7r + 7r _ )7r° 
events in the data (points). The hatched and open his- 
tograms represent the non-ISR background and sum of all 
backgrounds respectively. The smooth line approximates the 
non-ISR background, as described in the text. 

Figure 2 shows the 2(7r + 7r~)7r° invariant mass distribu- 
tion from threshold up to 4.5 GeV/c 2 for the experimental 
events in the signal region of Fig. 1. Narrow peaks are 
apparent at the J/tp and ip(2S) masses. The hatched his- 
togram represents the qq background, which is negligible 
at low mass but becomes large at higher masses. The 
open histogram represents the sum of all backgrounds, 
including those estimated from the control region. They 
total about 20% at low mass but account for 60-80% of 
the observed data near 4 GeV/c 2 . 

Considering uncertainties in the cross sections for the 
background processes, the normalization of events in the 
control region and the simulation statistics, we estimate 
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m(4rot°) (GeV/c 2 ) 

FIG. 3: (a) Invariant mass distributions for simulated 
2(7T + 7r~)7r° events in the phase space model, reconstructed 
in the signal (open) and control (hatched) regions of Fig. 1. 
(b) Net reconstruction and selection efficiency as a function 
of mass obtained from this simulation. The curve represents 
a 3 rd order polynomial fit. 



a systematic uncertainty on the signal yield that is about 
5% in the 1.5-1.8 GeV/c 2 mass region, but increases to 
20% at 2.5 GeV/c 2 and to more than 50% in the region 
above 3.5 GeV/c 2 . 



B. Selection Efficiency 

The selection procedures applied to the data are also 
applied to the simulated signal samples. The result- 
ing 2(7r + 7T~)7r° invariant-mass distributions in the signal 
and control regions of Fig. 1 are shown in Fig. 3(a) for 
the phase space simulation. The fraction of simulated 
events in the xl^o control region remains constant over 
mass, supporting the assumption of mass- independent x 2 
shape. The broad, smooth mass distribution is chosen to 
facilitate the estimation of the efficiency as a function 
of mass, and this model reproduces the observed distri- 
butions of pion momenta and polar angles. We divide 
the number of reconstructed simulated events in each 
25 MeV/c 2 mass interval by the number generated in 
that interval to obtain the efficiency shown as the points 
in Fig. 3(b); the curve represents a 3 rd order polynomial 
fit to the points. We simulate events with the ISR pho- 
ton confined to the angular range of EMC acceptance. 
The computed efficiency is for this fiducial region, but 
includes the acceptance for the final-state hadrons, the 
inefficiencies of the detector subsystems, and event loss 
due to additional soft-photon emission. 

The simulations including the LUir + Tr~ and/or r)Tr + Tr~ 
channels have very different distributions of mass and an- 



gles in the 2(ir + Tr~)ir° rest frame. However, the angular 
acceptance is quite uniform for ISR events, and the effi- 
ciencies are consistent with those from the phase space 
simulation within 3%. 




cm. 

FIG. 4: The e + e~— >2(7r + 7r~)7r () cross section as a function of 
e + e~ cm. energy measured with ISR data. Only statistical 
errors are shown. 

We study the shape of the X 2 W7r o distribution using 
events in the large J/ip peak. By comparing J/ip yields 
in data and simulation for x\ vn o 

< 40 and xl n7r o < 200, 
we limit any mis-modelling of the efficiency to 3%. We 
correct the track finding efficiency following the proce- 
dures described in Ref. [9], with a much larger sample 
of 2(7r + 7r~) events. We consider data and simulated 
events that contain a high-energy photon plus exactly 
three charged tracks and satisfy a set of kinematic crite- 
ria, including a good x 2 from a kinematic fit under the 
hypothesis that there is exactly one missing track in the 
event. We find that the simulated track-finding efficiency 
is overestimated by (0.8 ± 0.5)% per track, so we apply 
a correction of +(3 ± 2)% to the signal yield. We correct 
the 7r°-finding efficiency using the procedure described in 
detail in Ref. [10]. From ISR e + e~^cj7r°7^7r + 7r~7r 7r°7 
events selected with and without the 7r° from the u de- 
cay, we find an excess of simulated efficiency for one n° 
of (3 ± 2)%. 



C. Cross Section for e + e — > 2(7r + 7r )7T° 

We calculate the cross section as a function of effective 
cm. energy for the reaction e + e~—t 2(tt + -k~)tt® from 



Cr 2(TT+7T-)7T°{E C .m.) 



dN 2 

(7T+7T ) 7 r0 7 (£'c.m.) 

dC(E c _ m _) ■ e 2 ( 7r + 7r -) 7r o(.E c . m .) 



(2) 
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TABLE I: Measurements of the e + e — > 2(tt + tt )ty" cross section (errors are statistical only). 



£c.m. (GeV) 


a 


(nb) 


E e . m . (GeV) 


a 


(nb) 


Bcm. (GeV) 


a 


(nb) 


E e . m . (GeV) 


a 


(nb) 


1.0125 


0.02 


± 


0.03 


1.8875 


1.64 


± 


0.21 


2.7625 


0.39 


± 


0.13 


3.6375 


0.16 


± 0.08 


1.0375 


-0.01 


± 


0.04 


1.9125 


1.96 


± 


0.21 


2.7875 


0.20 


± 


0.12 


3.6625 


0.38 


± 0.09 


1.0625 


0.03 


± 


0.05 


1.9375 


1.88 


± 


0.21 


2.8125 


0.38 


± 


0.12 


3.6875 


1.55 


± 0.12 


1.0875 


-0.01 


± 


0.05 


1.9625 


1.76 


± 


0.20 


2.8375 


0.25 


± 


0.12 


3.7125 


0.40 


± 0.09 


1.1125 


0.04 


± 


0.06 


1.9875 


1.47 


± 


0.20 


2.8625 


0.07 


± 


0.12 


3.7375 


0.17 


± 0.08 


1.1375 


0.02 


± 


0.05 


2.0125 


1.66 


± 


0.20 


2.8875 


0.23 


± 


0.12 


3.7625 


0.14 


± 0.07 


1.1625 


0.09 


± 


0.07 


2.0375 


1.57 


± 


0.20 


2.9125 


0.30 


± 


0.12 


3.7875 


0.14 


± 0.07 


1.1875 


0.20 


± 


0.08 


2.0625 


1.77 


± 


0.20 


2.9375 


0.29 


± 


0.12 


3.8125 


0.08 


± 0.07 


1.2125 


0.13 


± 


0.09 


2.0875 


1.33 


± 


0.19 


2.9625 


0.35 


± 


0.12 


3.8375 


0.27 


± 0.07 


1.2375 


-0.02 


± 


0.10 


2.1125 


1.44 


± 


0.19 


2.9875 


0.45 


± 


0.11 


3.8625 


0.08 


± 0.06 


1.2625 


0.25 


± 


0.12 


2.1375 


1.42 


± 


0.19 


3.0125 


0.31 


± 


0.12 


3.8875 


0.09 


± 0.07 


1.2875 


0.31 


± 


0.13 


2.1625 


1.76 


± 


0.19 


3.0375 


0.39 


± 


0.12 


3.9125 


0.18 


± 0.07 


1.3125 


0.50 


± 


0.14 


2.1875 


1.38 


± 


0.18 


3.0625 


1.55 


± 


0.15 


3.9375 


0.15 


± 0.07 


1.3375 


0.72 


± 


0.16 


2.2125 


1.12 


± 


0.18 


3.0875 


17.13 ± 0.35 


3.9625 


0.14 


± 0.06 


1.3625 


0.91 


± 


0.17 


2.2375 


1.75 


± 


0.19 


3.1125 


10.76 ± 0.29 


3.9875 


0.03 


± 0.06 


1.3875 


0.95 


± 


0.18 


2.2625 


1.56 


± 


0.18 


3.1375 


0.90 


± 


0.13 


4.0125 


0.10 


± 0.06 


1.4125 


1.58 


± 


0.21 


2.2875 


1.20 


± 


0.17 


3.1625 


0.32 


± 


0.11 


4.0375 


0.11 


± 0.06 


1.4375 


1.65 


± 


0.22 


2.3125 


1.04 


± 


0.16 


3.1875 


0.23 


± 


0.11 


4.0625 


0.09 


± 0.06 


1.4625 


1.67 


± 


0.22 


2.3375 


1.25 


± 


0.17 


3.2125 


0.31 


± 


0.10 


4.0875 


0.05 


± 0.06 


1.4875 


1.92 


± 


0.24 


2.3625 


0.92 


± 


0.16 


3.2375 


0.25 


± 


0.10 


4.1125 


0.04 


± 0.06 


1.5125 


2.20 


± 


0.24 


2.3875 


0.82 


± 


0.15 


3.2625 


0.15 


± 


0.09 


4.1375 


0.13 


± 0.06 


1.5375 


2.20 


± 


0.24 


2.4125 


1.13 


± 


0.15 


3.2875 


0.17 


± 


0.10 


4.1625 


0.16 


± 0.06 


1.5625 


2.37 


± 


0.25 


2.4375 


0.58 


± 


0.14 


3.3125 


0.40 


± 


0.10 


4.1875 


0.18 


± 0.06 


1.5875 


2.36 


± 


0.25 


2.4625 


0.81 


± 


0.15 


3.3375 


0.07 


± 


0.09 


4.2125 


0.14 


± 0.06 


1.6125 


3.34 


± 


0.27 


2.4875 


0.64 


± 


0.15 


3.3625 


0.10 


± 


0.09 


4.2375 


0.11 


± 0.06 


1.6375 


3.29 


± 


0.27 


2.5125 


0.77 


± 


0.14 


3.3875 


0.24 


± 


0.09 


4.2625 


0.10 


± 0.06 


1.6625 


3.77 


± 


0.27 


2.5375 


0.52 


± 


0.14 


3.4125 


0.28 


± 


0.09 


4.2875 


0.05 


± 0.06 


1.6875 


3.20 


± 


0.26 


2.5625 


0.50 


± 


0.14 


3.4375 


0.15 


± 


0.09 


4.3125 


0.11 


± 0.05 


1.7125 


2.45 


± 


0.25 


2.5875 


0.42 


± 


0.13 


3.4625 


0.03 


± 


0.08 


4.3375 


0.08 


± 0.05 


1.7375 


2.71 


± 


0.24 


2.6125 


0.39 


± 


0.13 


3.4875 


0.13 


± 


0.08 


4.3625 


0.11 


± 0.05 


1.7625 


2.53 


± 


0.24 


2.6375 


0.65 


± 


0.14 


3.5125 


0.24 


± 


0.09 


4.3875 


0.08 


± 0.05 


1.7875 


2.31 


± 


0.23 


2.6625 


0.49 


± 


0.13 


3.5375 


0.06 


± 


0.08 


4.4125 


0.03 


± 0.05 


1.8125 


2.50 


± 


0.23 


2.6875 


0.51 


± 


0.13 


3.5625 


0.18 


± 


0.08 


4.4375 


0.11 


± 0.05 


1.8375 


2.14 


± 


0.22 


2.7125 


0.74 


± 


0.13 


3.5875 


0.16 


± 


0.08 


4.4625 


0.07 


± 0.05 


1.8625 


1.93 


± 


0.21 


2.7375 


0.29 


± 


0.12 


3.6125 


0.07 


± 


0.08 


4.4875 


0.15 


± 0.06 



ment of inclusive 2(7r + 7r~)7r° final state for a comparison. 
The applied corrections and systematic uncertainties are 
summarized in Tabic II. 



The cross section rises from threshold to a peak value 
of about 4.0 nb near 1.65 GeV, then generally decreases 
with increasing energy except for prominent peaks at 
the J/tp and ip(2S) masses. Gaussian fits to the sim- 
ulated line shapes give a resolution on the measured 
2(7r + 7T~)7r° mass that varies between 6.8 MeV/c 2 in 
the 1.5-2.5 GeV/c 2 region and 8.8 MeV/c 2 in the 2.5- 
3.5 GeV/c 2 region. The resolution function is not purely 
Gaussian due to soft-photon radiation, but less than 20% 
of the signal is outside the 25 MeV/c 2 mass bin. Since 
the cross section has no sharp structure other than the 
J/ip and ip(2S) peaks discussed in Sec. VIII below, we 
apply no correction for the resolution. 



where E c , m , = m 2 ( 7r + 7r -) 7I -oc 2 , m 2 ( 7r + 7r -) 7I -o is the measured 
invariant mass of the 2(7r + 7r~)7r° system, rf^ 2 ( 7r + T -) T o 7 
is the number of selected events after background sub- 
traction in the interval dE cm , and e2( 7r + 7r -) 7r o (E cm ) is 
the corrected detection efficiency. We calculate the dif- 
ferential luminosity, dC(E c , m ), in each interval dE c , m , 
using integrated BABAR luminosity and the probability 
density function from Eq. 1. We compare the experi- 
mental di-muon mass spectrum from ISR (i + [i~-f events, 
selected with the help of the instrumented flux return 
(IFR), with the calculated one [9] and conservatively esti- 
mate a systematic uncertainty on dC of 3%. This dC has 
been corrected for vacuum polarization (VP), so the ob- 
tained cross section includes contribution from VP which 
should be excluded when using these data in calculations 
of (g — 2) M [6]. The initial- and part of the final-state 
soft-photon emissions are canceled out in the ratio. 

We show the cross section as a function of energy in 
Fig. 4, with statistical errors only, and provide a list of 
our results in Table I. There is no direct e+e - measure- 
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FIG. 5: (a) The smallest tt + tt~tt mass in each selected FIG. 6: Scatter plots of (a) each tt + iy~ mass versus the mass 

2(7T + 7r~)7r° event versus the five-pion mass, (b) The full of the remaining 7r + 7r~7r° and (b) each tt^tt mass versus the 

7r + 7r _ 7T° mass distribution (four entries per event) in these mass of the remaining 7r + 7r~7r T for the selected 2(7r + 7r~)7r° 

events. The cross-hatched histogram represents the estimated candidates (four entries per plot per event). 
non-ISR background. 



D. Substructure in the 2(7r + 7r Final State 

The 2(7r + 7r _ )7T° final state has a rich internal structure. 
Figure 5(a) shows a scatter plot of the smallest 7r + 7r~7r° 
mass in each candidate event versus the five-pion mass. 
There are horizontal bands corresponding to the r]ir + Tr~ 
and W7r + 7r~ channels as well as vertical bands from the 
J /if) and i>(2S). Figure 5(b) shows the full tt + tt^tt° mass 
distribution (four entries per event) for selected events 
as the open histogram and for the estimated non-ISR 
background as the cross-hatched histogram. There is also 
a small signal for the 07r + 7r~ channel, and a peak at 
the J/ip mass, which is due to the tp(2S) — ► J/ipir + ir~ , 
J I 'tp — ► ir + tt~ ir° decay chain. 

Figure 6(a) shows a scatter plot of all four ■m^+ 1T - vs. 
m^+^-^o combinations in each event. There is a horizon- 
tal band corresponding to the p°(770) and an enhance- 
ment where it crosses the vertical 77 band. A / o ± (770) 
band is similarly visible in Fig. 6(b), a scatter plot of 
all four m 1T ± 7T o vs. m^+^-^ combinations in each event. 



TABLE II: Summary of systematic uncertainties on the 
e + e~ — ► 2(7T + 7r~)7r° cross section. The total uncertainty is 
the sum in quadrature of the components. 



Source Correction Uncertainty 



Rad. Corrections 




1% 








Backgrounds 




5% 


E c .m. 


=1.7 


GeV 






20% 


Ec.m. 


= 2.5 


GeV 






50% 


E c .m. 


= 3.5 


GeV 


Model Dependence 




3% 








xL^ Distn - 




3% 








Tracking Efficiency 


+3% 


2% 








7T° Efficiency 


+3% 


2% 








ISR Luminosity 




3% 








Total 


+6% 


7% 


E c .m. 


= 1.7 


GeV 






20% 


E c .m. 


= 2.5 


GeV 






50% 


E c .m. 


= 3.5 


GeV 



There is a suggestion of structure along these bands and 
in Fig. 5(b) around 1.2-1.3 GeV/c 2 , which could corre- 
spond to the oi(1260), 7r(1300) or ci2(1320) resonances. 
We now study events containing an 77, to or p in detail. 




m(jlVjl°) (GeV/c 2 ) 

FIG. 7: Distribution of the 7r + 7r~7r° mass closest to the r\ 
mass in the data (points). The histograms represent the dis- 
tributions from simulated r\p (open), tJ7r + 7r~ (dashed) and 
uds events (hatched), normalized as described in the text. 
The inner (inner and outer) vertical lines delimit the r]Tr + 7r~ 
signal (sideband) region. 



E. The r77r + 7T and r\p Intermediate States 

To extract the contribution from the ?y7r + 7r~ interme- 
diate state we select the 7r + 7r~7r° combination in each 
event (from four possible combinations) with mass clos- 
est to the 77 mass. Figure 7 shows the distribution of 
this mass in the data as points, along with various simu- 
lated distributions. The open histogram is for simulated 
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3 3.5 4 4.5 
m(r|7tV) (GeV/c 2 ) 

FIG. 8: Invariant mass distribution for events in the ?77r + 7r _ 
signal (open histogram) and sideband (dashed) regions (see 
Fig. 7). The points represent the background obtained from 
the xiLir control region, and the hatched histogram is that 
from uds events. The inset shows the low end of the mass 
distribution, where direct e + e~— *r} lr ) events are visible. 




m(7tV) (GeV/c 2 ) 

FIG. 9: Invariant mass distribution of the 7r + 7r _ pair not 
from the r\ in selected t]tt + it~ events in the data (points) and 
in simulated r\p events (histogram). 



i]p— >2(7r + 7r~)7r° events, is normalized to data, and shows 
only a narrow 77 peak. The dashed histogram for simu- 
lated LJ7r + 7r~ events shows a strong uj peak with a tail 
toward lower masses that contributes a small number of 
events in the r\ region. The hatched histogram for simu- 
lated uds events is normalized as described in Sec. IV A 



and shows both 77 and u signals over a small combinatoric 
contribution. 

We define an rj signal region as mass in the range 525- 
575 MeV/c 2 , indicated by the inner vertical lines in Fig. 7, 
and two sidebands, 500-525 and 575-600 MeV/c 2 , indi- 
cated by the outer vertical lines. The 77 signal region con- 
tains 1897 data events, and we show their 2(7r + 7r _ )7r° in- 
variant mass distribution as the open histogram in Fig. 8. 
The hatched histogram is the uds background, evaluated 
as described in section IV A, which contributes mostly at 
higher masses. We evaluate the remaining background 
in two ways. Using the control region of the xt^n 
tribution, as discussed in Sec. IV A but only in the rj 
signal region, we obtain the mass distribution shown as 
the points in Fig. 8. Alternatively, the mass distribution 
for events in the rj sidebands is shown as the dashed his- 
togram. These two distributions are consistent, indicat- 
ing that very few non-77 2(7r + 7r~)7r° events are present. 
Here we use the sideband distribution since it is more 
precise and contains all backgrounds. The inset in Fig. 8 
shows the distribution in the region below 1.3 GeV/c 2 
with finer binning. A signal from the direct e + e~— >rj'^/ 
process is visible; these events were studied in our previ- 
ous measurement [25]. 

The invariant mass distribution of the 7r + 7r~ pair not 
from the 77 is shown after subtraction of the uds and 77- 
sideband backgrounds as the points in Fig. 9, and has a 
strong peak in the p(770) region. The histogram in Fig. 9 
is the distribution for simulated rjp — > 2(7r + 7r _ )7r° events, 
and its similarity to the data indicates that this channel 
dominates the 7i7r + 7r~ intermediate state. We therefore 
use the simulated r\p events to estimate the detection 
efficiency for r]ir + ir~ events, although the other simu- 
lations give consistent results. Figure 10(a) shows the 
simulated invariant mass distribution for selected events, 
and Fig. 10(b) shows the simulated efficiency, which in- 
cludes the 22.6% branching fraction of 77 — >7r + 7T~7r°. 

Subtracting the backgrounds and dividing by the ISR 
luminosity and efficiency, parametrized by the third or- 
der polynomial fit shown in Fig. 10(b) and corrected as 
discussed above, we obtain the e + e~^?77r + 7r~ cross sec- 
tion shown in Fig. 11. Also shown are the previous direct 
e+e" measurements from the DM2 [13], CMD2 [15] and 
ND [16] experiments. All measurements are consistent, 
and ours covers the widest energy range and is by far the 
most precise above 1.4 GeV. 

The cross section shows a steep rise from r]p(770) 
threshold, followed by a general decrease with increasing 
energy. Possible structures near 1.6 and 1.8 GeV can- 
not be resolved with the current statistics. We list the 
cross section in Table III for cm. energies up to 3 GeV 
with statistical errors only. The systematic uncertainties 
are the same as those discussed in section IV C, totalling 
about 8% below 3 GeV. Above 3 GeV the cross section is 
consistent with zero within the current statistical errors, 
except for the J/ip peak, which is discussed below. 
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TABLE III: Measurements of the e + e — » r/7r + 7r cross section (errors are statistical only). 



E c . m . (GeV) 


a (nb) 


£c.m. (GeV) 


a (nb) 


E c .m. (GeV) 


a (nb) 


E c .m. (GeV) 


a (nb) 


1.0250 


0.00 ± 0.05 


1.5250 


4.29 ± 0.37 


2.0250 


0.54 ± 0.14 


2.5250 


0.00 ± 0.07 


1.0750 


0.11 ± 0.08 


1.5750 


3.13 ± 0.34 


2.0750 


0.50 ± 0.15 


2.5750 


0.06 ± 0.07 


1.1250 


0.05 ± 0.10 


1.6250 


2.83 ± 0.35 


2.1250 


0.36 ± 0.12 


2.6250 


0.11 ± 0.07 


1.1750 


0.20 ± 0.13 


1.6750 


2.13 ± 0.29 


2.1750 


0.64 ± 0.12 


2.6750 


0.02 ± 0.07 


1.2250 


0.10 ± 0.13 


1.7250 


2.33 ± 0.28 


2.2250 


0.32 ± 0.11 


2.7250 


0.05 ± 0.06 


1.2750 


0.66 ± 0.18 


1.7750 


1.90 ± 0.25 


2.2750 


0.42 ± 0.10 


2.7750 


0.00 ± 0.06 


1.3250 


0.76 ± 0.24 


1.8250 


1.57 ± 0.23 


2.3250 


0.16 ± 0.10 


2.8250 


0.05 ± 0.06 


1.3750 


1.37 ± 0.26 


1.8750 


0.80 ± 0.18 


2.3750 


0.15 ± 0.09 


2.8750 


0.00 ± 0.05 


1.4250 


2.88 ± 0.33 


1.9250 


0.77 ± 0.17 


2.4250 


0.05 ± 0.08 


2.9250 


0.20 ± 0.06 


1.4750 


3.59 ± 0.34 


1.9750 


0.69 ± 0.15 


2.4750 


0.25 ± 0.09 


2.9750 


0.00 ± 0.05 
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FIG. 10: (a) Invariant mass distribution for selected simulated 
rjp events, and (b) the detection efficiency vs. mass, including 
the rj^ 7r + 7T-~7r (J branching fraction. 



+ 
+ 




FIG. 11: The e e~ — > rjn ty~ cross section as a function of 
cm. energy obtained via ISR at BABAR. The direct mea- 
surements from DM2, ND and CMD2 are also shown. Only 
statistical errors are shown. 



F. The lj7t 7T and ujfo Intermediate States 

To extract the contribution of the c<j7r + 7r~ intermedi- 
ate state we select the 7r + 7r~7r° combination with mass 
closest to the ui mass. Events with this mass below 
0.6 GeV/c 2 are predominantly r/ir + ir~ events, and we ig- 
nore them. Subtracting the simulated uds background 
and ISR-type background from the X 2 ^ control region, 
as described in Sec. IV A, we obtain the mass distribution 
shown in Fig. 12 as points. The histogram is for simu- 
lated oj7r + 7r~ events and describes the peak in the data 
well, but a background from non-w7r + 7r~, non-?77r + 7r _ 
events is still present. 

We define an uj signal region as 7r + 7r~7r° mass in the 
range 745-825 MeV/c 2 , indicated by the inner vertical 
lines in Fig. 12 and containing 7693 events, and two side- 
bands, 706-745 and 825-865 MeV/c 2 , indicated by the 
outer vertical lines. Figure 13 shows the invariant mass 



distributions for events in the signal region (points) and 
sidebands (hatched histogram). The sideband events are 
ISR 2(7r + 7T~)7r° events but without an u> or r/. They 
contribute mostly at higher energies including the J /tp 
peak. 

We evaluate the detection efficiency using the (jjtt + ti~ 
phase space simulation. It is similar to that in Fig. 3, 
differing by few percent due to the additional selection 
criteria. Subtracting the sideband background and di- 
viding by the corrected efficiency, ISR luminosity and 
the 89.1% branching fraction of u — > ir + Tr~ir°, we ob- 
tain the e + e~— > lutt + tt~ cross section shown in Fig. 14. 
Also shown are the previous direct e + e~ measurements 
from the DM2 [14], DM1 [12] and CMD2 [15] experi- 
ments. All measurements are consistent, and ours cover 
the widest energy range and are by far the most precise 
above 1.4 GeV. 
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TABLE IV: Measurement of the e+e" 



cross section (errors are statistical only). 



E c . m . (GeV) 


a 


(nb) 


£c.m. (GeV) 


a 


(nb) 


£c.m. (GeV) 


a 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


1.1500 


0.00 


± 


0.09 


1.5000 


1.06 


± 


0.25 


1.8500 


0.79 


± 


0.20 


2.2000 


0.11 


± 0.14 


1.1750 


0.00 


± 


0.05 


1.5250 


1.33 


± 


0.24 


1.8750 


0.84 


± 


0.20 


2.2250 


0.52 


± 0.14 


1.2000 


0.06 


± 


0.08 


1.5500 


1.67 


± 


0.27 


1.9000 


0.95 


± 


0.20 


2.2500 


0.35 


± 0.13 


1.2250 


0.00 


± 


0.12 


1.5750 


1.30 


± 


0.27 


1.9250 


0.61 


± 


0.19 


2.2750 


0.27 


± 0.12 


1.2500 


0.15 


± 


0.14 


1.6000 


2.10 


± 


0.28 


1.9500 


0.28 


± 


0.18 


2.3000 


-0.04 


± 0.12 


1.2750 


0.23 


± 


0.14 


1.6250 


2.21 


± 


0.29 


1.9750 


0.69 


± 


0.17 


2.3250 


0.11 


± 0.11 


1.3000 


0.30 


± 


0.15 


1.6500 


2.80 


± 


0.30 


2.0000 


0.11 


± 


0.17 


2.3500 


-0.03 


± 0.10 


1.3250 


0.33 


± 


0.16 


1.6750 


2.19 


± 


0.28 


2.0250 


0.49 


± 


0.16 


2.3750 


0.00 


± 0.10 


1.3500 


0.55 


± 


0.20 


1.7000 


1.99 


± 


0.26 


2.0500 


0.70 


± 


0.15 


2.4000 


0.05 


± 0.10 


1.3750 


0.88 


± 


0.20 


1.7250 


1.38 


± 


0.25 


2.0750 


0.40 


± 


0.15 


2.4250 


0.02 


± 0.10 


1.4000 


0.69 


± 


0.22 


1.7500 


1.51 


± 


0.24 


2.1000 


0.35 


± 


0.15 


2.4500 


-0.04 


± 0.09 


1.4250 


0.83 


± 


0.24 


1.7750 


1.45 


± 


0.23 


2.1250 


0.26 


± 


0.14 


2.4750 


0.05 


± 0.10 


1.4500 


1.17 


± 


0.23 


1.8000 


1.18 


± 


0.23 


2.1500 


0.56 


± 


0.14 


2.5000 


0.02 


± 0.10 


1.4750 


0.95 


± 


0.25 


1.8250 


1.19 


± 


0.21 


2.1750 


0.43 


± 


0.14 


2.5250 


0.15 


± 0.09 




0.8 0.9 1 

m(% + %%°) (GeV/c 2 ) 

FIG. 12: Distribution of the 7r + 7r - 7r° mass closest to the w 
mass in the data (points) after subtraction of uds and ISR- 
type backgrounds. The histogram is the distribution for sim- 
ulated u)it + -k~ events. The inner (inner and outer) vertical 
lines delimit the uj signal (sideband) region. 



The cross section is consistent with zero below 
1.2 GeV, then rises to a peak value of about 2.5 nb at 
about 1.65 GeV, followed by a general decrease with in- 
creasing energy. We list the cross section in Table IV for 
cm. energies up to 2.4 GeV with statistical errors only. 
The systematic uncertainties are the same as those dis- 
cussed in section IV C, totalling about 8% in this range. 
Above 2.4 GeV the cross section is consistent with zero 
within the current statistical errors, except for the J/tp 
peak, which is discussed below. 

For events in the w signal region with a five-pion mass 
below 3.0 GeV/c 2 , we show the invariant mass distri- 




1.5 2 2.5 3 3.5 4 ,4.5 
m(awtV) (GeV/c 2 ) 

FIG. 13: Invariant mass distribution for events in the lu-k + it~ 
signal region (points) and sidebands (hatched histogram). 



bution of the 7r + 7r~ pair not from the u) in Fig. 15(a). 
A peak is visible in the data at the Jo (980) mass, and 
the histogram is for a simulation that includes ujk + -k~ 
phase space and w/o(980) combined so as to describe the 
data. We define an w/o(980) signal region by this 7r + 7r~ 
mass in the region 0.88-1.04 GeV/c 2 , indicated by the in- 
ner vertical lines in Fig. 15(a), and sidebands 0.80-0.88 
and 1.04-1.12 GeV/c 2 , indicated by the outer lines. Sub- 
tracting the sideband contribution from that in the signal 
region and dividing by the corrected efficiency, ISR lumi- 
nosity and the 2/3 branching fraction of /o(980)^7r + 7r~ 
(assuming 27r decay mode dominance [5]), we obtain the 
e + e~ — > cj/ (980) cross section shown in Fig. 15(b) and 
listed in Table V. This measurement of the cross sec- 
tion shows a very fast rise from threshold and a possible 
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TABLE V: Measurements of the e + e — » lo /o(980) cross section (errors are statistical only). 



E c . m . (GeV) 


a (nb) 


E c . m . (GeV) 


a (nb) 


E c . m . (GeV) 


a (nb) 


E c .m. (GeV) 


a (nb) 


1.7000 


0.02 ± 0.14 


1.9000 


0.51 ± 0.18 


2.1000 


0.39 ± 0.11 


2.3000 


0.06 ± 0.07 


1.7250 


0.02 ± 0.14 


1.9250 


0.48 ± 0.16 


2.1250 


0.18 ± 0.11 


2.3250 


0.15 ± 0.06 


1.7500 


0.10 ± 0.18 


1.9500 


0.32 ± 0.15 


2.1500 


0.31 ± 0.11 


2.3500 


0.02 ± 0.05 


1.7750 


0.79 ± 0.17 


1.9750 


0.42 ± 0.15 


2.1750 


0.14 ± 0.10 


2.3750 


0.03 ± 0.05 


1.8000 


0.89 ± 0.20 


2.0000 


0.17 ± 0.12 


2.2000 


0.17 ± 0.09 


2.4000 


0.06 ± 0.05 


1.8250 


0.91 ± 0.17 


2.0250 


0.34 ± 0.12 


2.2250 


0.11 ± 0.11 


2.4250 


0.07 ± 0.05 


1.8500 


0.45 ± 0.16 


2.0500 


0.27 ± 0.10 


2.2500 


0.03 ± 0.09 


2.4500 


0.01 ± 0.04 


1.8750 


0.46 ± 0.16 


2.0750 


0.27 ± 0.12 


2.2750 


0.02 ± 0.08 


2.4750 


0.01 ± 0.05 
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E c . m . (GeV) 

FIG. 14: The e + e~ — > lott + tt~ cross section as a function of 
cm. energy obtained via ISR at BABAR,. The direct measure- 
ments from DM1, DM2 and CMD2 are also shown. Only 
statistical errors are shown. 



structure at about 1.85 GeV/c 2 , followed by a monotonic 
decrease with increasing energy. 

We subtract the e + e~ — > ujf {980) cross section from 



the inclusive e + e 



cross section to obtain the 



cross section shown in Fig. 16. A peak is visible, presum- 
ably from the w(1650), and the shoulder at lower masses 
can be attributed to the w(1420). We fit this cross sec- 
tion as a function of E c m = ^fs from threshold up to 
2.4 GeV with a sum of vector resonances, 



P(s) 



^okmlTk 



fe=i 



il-s- i^fsT k y/Ptmgj 



12-KBkeBkfC 



(3) 
(4) 



P{s) = /s 



where Tk and rrik are the full width and mass of the 
fc th resonance, Bj~ e and B k f are its branching fractions 




0.5 I 



m(;iV) (GeV/c*) 



(GeV) 



FIG. 15: (a) Invariant mass distribution of the tt + tt~ pair not 
from the lo in selected lutv + tv~ events in the data (points) and 
in simulated events (histogram). The vertical lines delimit the 
/o(980) signal region and sidebands, (b) The e + e~— ^w/o(980) 
cross section. 



to e + e and the final state / 



respectively, 



P(s) is a simple approximation of the phase space with 
a threshold cutoff at m = 1.2 GeV, and C — 3.893 • 
10 5 nb GeV 2 is a conversion constant. This formulation 
allows the extraction not only of the product B ke B k f but 
also the peak cross section a ok or the product T ke B v f. 

We consider the resonances, k = ui(782) (lj), w(1420) 
(lo 1 ), and w(1650) (to"), where w(782) is below thresh- 
old and is used as a convenient "coherent background" 
with parameters fixed to PDG values [5] and with </> w (782) 
set to 0. We perform three fits, the results of which are 
listed in Table VI and compared with results of a similar 
fit from our study of the ISR 7r + 7r~7r° process [8] and 
with current PDG values [5]. In the first fit, we set the 
contribution from uj to zero and set = 0. The fit- 
ted cross section is dominated by lo" , and the lo' has a 
relatively narrow width. 

Next we float the contribution from lo but fix the 
relative phases to the values used in our ISR 7r + 7r~7r° 
study [8], (j)^' = 7r and <fi u " =0. The resulting contribu- 
tion from lo' is almost 10 times higher due to destruc- 
tive interference, but the masses and widths are similar 
to those from the first fit. In particular the lo' width is 
lower than that found in Ref . [8] . The fitted peak lo cross 
section corresponds to a large a>7r + 7r~ branching fraction 
of about 7%, but this is driven by the data above 2 GeV 
and should be considered a measure of the coherent back- 
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TABLE VI: Summary of parameters obtained from the fits described in the text. The values without errors were fixed in that 
fit. 



Fit 


1 


2 


3 


3-7T [8] 


PDG [5] 


a 0w , (nb) 


0.10±0.08 


1.01±0.29 


0.64±0.34 






B ee B w / f ■ 10 6 


0.013±0.010 


0.13±0.04 


0.101±0.055 


0.82±0.08 


— 


T ee B w , f (eV) 
m w ,(GeV/c 2 ) 


1.4±1.01 


17.5±5.4 


37.8±12.1 


369 


— 


1.381±0.032 


1.382±0.023 


1.463±0.070 


I.350±0.030 


1.40-1.45 


T w > ( GeV) 


0.105±0.090 


0.133±0.048 


0.383±0.233 


0.450±0.140 


0.180-0.250 


<j> w i (rad.) 


-1.93±0.73 


TT 


-0.61±0.94 


TT 




a 0w ii (nb) 


2.14±0.18 


2.47±0.18 


1.03±0.54 






B ee B w rif ■ 10 6 


0.41±0.03 


0.47±0.04 


0.193±0.087 


1.3±0.2 




r ee B w „ f (eV) 


96.5±10.9 


103.5±8.3 


28.7±7.7 


286 




m w ii ( GeV/c ) 


1.673±0.011 


1.667±0.013 


1.661±0.032 


1.660±0.011 


1.670±0.030 


IV (GeV) 


0.236±0.029 


0.222±0.025 


0.148±0.037 


0.220±0.040 


0.315±0.035 


<j> w ii (rad.) 








0.02±0.71 







a 0w (nb) 





102±67 


I47±140 


PDG 




X7n-d.f. 


36.2/48 


34.9/48 


32.2/46 
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FIG. 16: The e + e~ — > lutt + tt~ cross section excluding the 
u/o(980) contribution (points). The curve shows the result 
of the fit of the w(1420) and w(1650) resonances described in 
the text (Fit 3 in Table VI). 



ground. 

An interference with other unaccounted vector mesons 
could produce deviations from the assumed values of the 
phases. To demonstrate this we float the uu level and 
both relative phases, and show the result as the curve in 
Fig. 16. The coherent background is larger and both the 
lo' and lo" peak cross sections are much lower than in the 
second fit. Both masses are consistent with the other fits, 
but the lo 1 {lo") width is larger (smaller) and statistically 
consistent with our ISR tt + tt~tt° study. 

A better understanding of the background, including 
any structure above 2 GeV and any contribution from 
excited p or <j> states, is needed in order to make pre- 



cise measurements of the excited lo resonance parame- 
ters. Taking the results from the second fit and using 
the differences from the other fits to estimate systematic 
errors, we obtain: 



"^(1420) 
r w (1420) 

^w(1650) 
^(1650) 



1.38 ±0.02 ±0.07 GeV/c 2 , 
0.13 ±0.05 ±0.10 GeV, 
1.667 ±0.013 ±0.006 GeV/c 2 , 
0.222 ±0.025 ±0.020 GeV. 



The w(1650) width is significantly different from the PDG 
value [5] based on the DM2 results [14, 15], but consistent 
with our measurement in ISR 7r + 7r _ 7r° events [8]. Note 
that the structure, observed in our study of ISR w(782)?7 
events [10] and described by a resonance with m = 1.645± 
0.008 GeV/c 2 and T = 0.114 ± 0.014 GeV can also be 
interpreted as w(1650). 



G. The ,9(770)377 Intermediate States 

To study events containing a charged or neutral p(770) 
we first exclude any event in which a 7r + 7r~7r° combina- 
tion has invariant mass within 25 MeV/ c 2 of the r\ mass 
or within 40 MeV/ c 2 of the lo mass. For this study we also 
exclude events with a five-pion mass within 50 McV/c 2 of 
the J/ip mass. Figure 17 shows the invariant mass distri- 
butions for all four tt + tt~ pairs and all four tt^tt pairs 
in the remaining events. The ISR and non-ISR back- 
grounds are subtracted using the procedures described 
above. These two distributions are quite similar and 
show strong p(770) peaks. The hatched histogram in 
Fig. 17 shows the mass distribution for the n + ir + and 
7r~7r~ pairs (two entries per event), which gives an es- 
timate of the combinatorial background. The difference 
between these distributions is consistent with an average 
of two p per event: one p is charged and the other neu- 
tral, since the yields are consistent and e + e~^> p° p a n° is 
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FIG. 17: The background-subtracted invariant mass distri- 
butions for all n + n~ pairs (open histogram, four entries per 
event), 7r ± 7r° pairs (points, four entries per event) and ty + ty + 
and 7r~7r~ pairs (hatched histogram, two entries per event), 
for events with no r\ or oj-^7r + 7r~7r° candidates. 




m(37t) (GeV/c 2 ) ! m(3it) (GeV/c 2 ) 



FIG. 18: (a) The invariant mass distributions for 7r + 7r~7r 
(triangles) and 7r + 7r~7r ± (squares) combinations for events in 
which the other 7r + 7r~ or -k^h pair is in the p signal region. 
The hatched histogram is an estimate of the combinatorial 
background, (b) The average of the triangles and squares 
in (a) minus the combinatorial background. The line is the 
result of a Breit-Wigner fit. 



forbidden by C-parity. This suggests one or more quasi- 
two-body intermediate states, X ±,0 p T '°, where X could 
be ai(1260), tt(1300) or a 2 (1370), which have I = 1 and 
a dominant pn decay. 

We now select events that contain a 7r ± 7r° or tt + tt~ 
pair with mass within 150 MeV/c 2 of the p mass. Fig- 
ure 18(a) shows the mass distributions for the other 
charged (squares) and neutral (triangles) three-pion com- 
binations (up to four total entries per event). These two 
distributions are consistent and the hatched histogram is 
an estimate of the combinatorial background from dou- 
bly charged (7T + 7r + 7r° and 7r~7r~7r°) combinations. By 
averaging the charged and neutral distributions and sub- 




FIG. 19: (a) The total e + e" ^2(-k + ty~)-k° cross section (cir- 
cles) and the contributions from lutv + iv~ (squares) and r]Tr + n~ 
(triangles), (b) The cross section obtained as the difference 
between the total and the latter two contributions, which is 
dominated by the e + e~— +p(770)X process. 



trading the combinatorial background we obtain the dis- 
tribution shown in Fig. 18(b), which is consistent with a 
resonant structure. Fitting a single Breit-Wigner func- 
tion gives 

m(X) = 1.243 ±0.012 ±0.020 GeV/c 2 ; 
T(X) = 0.410 ±0.031 ±0.030 GeV. 

The first errors are statistical and the second system- 
atic, dominated by the background subtraction proce- 
dure. These values are inconsistent with the 02 (1320) res- 
onance, but consistent, within large uncertainties, with 
the 7r(1300) and ai(1260) [5]. An angular analysis could 
distinguish between these possibilities, but requires sub- 
stantially higher statistics due to the large combinatorial 
background. 

Since the events that do not contain an 77 or oj appear 
to be predominantly Xp(77Q) — > p°p ± iT T events, where 
X is consistent with a single resonance, we obtain an 
e+e~—> Vp(770) cross section as the difference between 
the total e + e~^ 2(7r + 7r~)7r° cross section (Fig. 4) and 
the e + e~ — > W7r + 7r~ and e + e~ — > r)ir + Tr~ cross sections 
(Figs. 14 and 10 with branching fraction corrections re- 
moved). We show these three cross sections in Fig. 19(a) 
for energies up to 3 GeV, and the difference in Fig. 19(b); 
it shows no sharp structure. Above 3 GeV, the contribu- 
tions from 777r + 7r~ and lutt + tt~ are consistent with zero, 
so the cross section is as in Fig. 4, except for the J/tp and 
ip{2S) peaks, which can have a different substructure (see 
Sec. VIII). 

V. THE 2(7T+7r-)r7 FINAL STATE 
A. Final selection and background 

To suppress K + K~ n + ix~ r\ background, we require 
that no more than one track in the event is identi- 
fied as a kaon, and we also fit under all four possi- 
ble K + K~n + n~ri hypotheses and require X^a™ > ^0- 
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Signal 



Control 




x (4jcn) 



FIG. 20: Distribution of \ from the 5C fit for 2(7r + 7r")?7 
candidates in the data (points). The open and hatched his- 
tograms are the distribution for simulated signal events and 
non-ISR background, respectively, normalized as described in 
the text. 



We suppress any K ± K^-k^t] background by requiring 
all tracks to extrapolate within 2.5 mm of the beam 
axis. The xl^?? distribution for the remaining events is 
shown as points in Fig. 20, and the distribution for simu- 
lated 2(7r+7r- ~)ir° events (open histogram) is normalized 
to the data in the region xt,™ < 10- We do not simulate 
2(7r + 7r~)?7 events, but we expect the efficiency and res- 
olution for r\ to be similar to that for 7r°. The hatched 
histogram represents the non-ISR background contribu- 
tion obtained from the JETSET simulation. It is domi- 
nated by 2(7r + 7r~)7r°?7 events, and we use the same nor- 
malization factor as for the 2(7r + 7r _ )7r° events described 
in Sec. IV B. We define a signal region, x| < 40, con- 
taining 4272 events, and a control region for the estima- 
tion of other backgrounds, 40 < xl^r; < 80, containing 
1485 events. 




m(2(7tV)Ti) (GeV/c ) 



FIG. 21: Invariant mass distribution for selected 2(7r + 7r - )?7 
events in the data (points). The hatched and open his- 
tograms represent, cumulatively, the non-ISR background and 
the background from the control region of Fig. 20. 




4 

E (GeV) 

cm. v ' 

FIG. 22: The e + e~— > 2(7r + 7r~)j7 cross section as a function of 
cm. energy measured with ISR data. Only statistical errors 
are shown. 



Figure 21 shows the 2(7r + 7r _ )7y invariant mass distri- 
bution from threshold up to 4.5 GeV/c 2 for events in 
the signal region. A J/ip signal is visible. The hatched 
histogram represents the non-ISR background, and the 
open histogram represents the sum of all backgrounds, 
where the ISR-type background is estimated from the 
control region. Both backgrounds are relatively small at 
low mass, about 20% altogether, but they account for 
50-80% of the observed data in the 3.0-4.5 GeV/c 2 re- 
gion. We subtract this total background in each bin to 
obtain a number of signal events. 



B. Cross section for e + e — > 2(7T + 7r )r] 

We calculate the cross section for the e + e~— > 2(ir + ir~ )rj 
process as described in Sec. IV C, by dividing the number 
of events in each 2(ir + ir~)ri mass bin by the corrected de- 
tection efficiency and differential luminosity. The angular 
acceptance is uniform in all of our simulations, and this 
has been demonstrated in our previous studies of four- 
and six-pion final states [9, 10]. We therefore use the 
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TABLE VII: Measurements of the e + e — > 2{tx + -k )r] cross section (errors are statistical only). 



m . (GeV) 


a 


(nb) 


E c .m. (GeV) 


a 


(nb) 


E c . m . (GeV) 


a 


(nb) 


£c.m. (GeV) 


a (nb) 


1.3125 


0.00 


± 


0.00 


2.1125 


0.84 


± 


0.18 


2.9125 


0.30 


± 


0.13 


3.7125 


0.21 ± 0.08 


1.3375 


0.04 


± 


0.04 


2.1375 


1.10 


± 


0.18 


2.9375 


0.52 


± 


0.14 


3.7375 


0.10 ± 0.07 


1.3625 


0.00 


± 


0.00 


2.1625 


1.03 


± 


0.20 


2.9625 


0.33 


± 


0.14 


3.7625 


0.04 ± 0.08 


1.3875 


0.04 


± 


0.04 


2.1875 


1.23 


± 


0.19 


2.9875 


0.69 


± 


0.14 


3.7875 


0.02 ± 0.08 


1.4125 


0.11 


± 


0.07 


2.2125 


0.96 


± 


0.18 


3.0125 


0.13 


± 


0.13 


3.8125 


0.13 ± 0.08 


1.4375 


0.05 


± 


0.06 


2.2375 


1.15 


± 


0.19 


3.0375 


0.65 


± 


0.14 


3.8375 


0.06 ± 0.07 


1.4625 


0.00 


± 


0.07 


2.2625 


1.21 


± 


0.18 


3.0625 


0.43 


± 


0.13 


3.8625 


0.13 ± 0.08 


1.4875 


0.08 


± 


0.07 


2.2875 


0.85 


± 


0.17 


3.0875 


1.45 


± 


0.19 


3.8875 


0.06 ± 0.07 


1.5125 


0.07 


± 


0.05 


2.3125 


1.00 


± 


0.17 


3.1125 


0.89 


± 


0.16 


3.9125 


0.18 ± 0.08 


1.5375 


0.01 


± 


0.07 


2.3375 


1.25 


± 


0.18 


3.1375 


0.47 


± 


0.12 


3.9375 


0.02 ± 0.07 


1.5625 


-0.06 


± 


0.06 


2.3625 


1.28 


± 


0.18 


3.1625 


0.48 


± 


0.12 


3.9625 


0.15 ± 0.08 


1.5875 


0.06 


± 


0.07 


2.3875 


0.74 


± 


0.17 


3.1875 


0.48 


± 


0.14 


3.9875 


0.10 ± 0.06 


1.6125 


-0.01 


± 


0.08 


2.4125 


0.99 


± 


0.17 


3.2125 


0.58 


± 


0.12 


4.0125 


0.10 ± 0.07 


1.6375 


-0.08 


± 


0.08 


2.4375 


0.92 


± 


0.17 


3.2375 


0.39 


± 


0.12 


4.0375 


-0.02 ± 0.06 


1.6625 


0.08 


± 


0.09 


2.4625 


0.91 


± 


0.17 


3.2625 


0.25 


± 


0.11 


4.0625 


0.11 ± 0.06 


1.6875 


0.09 


± 


0.08 


2.4875 


0.95 


± 


0.17 


3.2875 


0.31 


± 


0.11 


4.0875 


0.10 ± 0.06 


1.7125 


0.18 


± 


0.10 


2.5125 


0.68 


± 


0.16 


3.3125 


0.16 


± 


0.09 


4.1125 


-0.02 ± 0.05 


1.7375 


-0.04 


± 


0.08 


2.5375 


0.80 


± 


0.17 


3.3375 


0.12 


± 


0.11 


4.1375 


0.04 ± 0.07 


1.7625 


0.22 


± 


0.12 


2.5625 


0.86 


± 


0.16 


3.3625 


0.14 


± 


0.09 


4.1625 


0.05 ± 0.07 


1.7875 


0.31 


± 


0.13 


2.5875 


0.63 


± 


0.16 


3.3875 


0.19 


± 


0.11 


4.1875 


0.06 ± 0.05 


1.8125 


0.39 


± 


0.11 


2.6125 


0.53 


± 


0.15 


3.4125 


0.12 


± 


0.10 


4.2125 


0.06 ± 0.07 


1.8375 


0.63 


± 


0.16 


2.6375 


0.65 


± 


0.14 


3.4375 


0.42 


± 


0.11 


4.2375 


0.16 ± 0.06 


1.8625 


0.31 


± 


0.13 


2.6625 


0.81 


± 


0.15 


3.4625 


0.32 


± 


0.10 


4.2625 


0.04 ± 0.06 


1.8875 


0.27 


± 


0.14 


2.6875 


0.73 


± 


0.16 


3.4875 


0.13 


± 


0.10 


4.2875 


0.03 ± 0.07 


1.9125 


0.69 


± 


0.15 


2.7125 


0.43 


± 


0.13 


3.5125 


0.12 


± 


0.09 


4.3125 


0.04 ± 0.06 


1.9375 


0.59 


± 


0.17 


2.7375 


0.46 


± 


0.14 


3.5375 


0.19 


± 


0.09 


4.3375 


0.09 ± 0.06 


1.9625 


0.50 


± 


0.15 


2.7625 


0.64 


± 


0.15 


3.5625 


0.22 


± 


0.09 


4.3625 


-0.01 ± 0.07 


1.9875 


0.88 


± 


0.17 


2.7875 


0.65 


± 


0.15 


3.5875 


0.22 


± 


0.09 


4.3875 


-0.06 ± 0.06 


2.0125 


0.54 


± 


0.17 


2.8125 


0.61 


± 


0.14 


3.6125 


0.13 


± 


0.08 


4.4125 


-0.03 ± 0.06 


2.0375 


0.87 


± 


0.19 


2.8375 


0.39 


± 


0.13 


3.6375 


0.12 


± 


0.08 


4.4375 


0.10 ± 0.06 


2.0625 


1.32 


± 


0.20 


2.8625 


0.54 


± 


0.14 


3.6625 


0.19 


± 


0.08 


4.4625 


0.13 ± 0.07 


2.0875 


0.99 


± 


0.19 


2.8875 


0.65 


± 


0.14 


3.6875 


0.34 


± 


0.09 


4.4875 


0.07 ± 0.06 



same detection efficiency as for the 2(7r + 7r _ )7r° process, 
shown in Fig. 3, divided by the 77—^77 branching fraction 
of 39.28% [5], and with the systematic error increased to 
5%. We use the same corrections and uncertainties for 
the x 2 cut, tracking efficiency and rj- finding efficiency. 

We show the cross section as a function of energy in 
Fig. 22 with statistical errors only, and provide a list of 
our results in Table VII. This is the first measurement 
of this cross section, which shows a peak value of about 
1.2 nb at about 2.2 GeV, followed by a monotonic de- 
crease toward higher energies, broken only by a peak at 
the J/i/j mass, discussed in Sec. VIII. Again, the energy 
resolution is much smaller than the bin width and we ap- 
ply no correction. The overall systematic error is about 
10% for energies below 3 GeV, rising to 30-50% in the 
3-4.5 GeV region. 

C. Substructure in the 2(7r~'~7r _ )r) Final State 

We might expect a rich internal structure in the 
2(7r + 7r - )?7 final state. Figure 23(a) shows the invari- 
ant mass distribution of the 7r + 7r~7r + 7r~ system recoiling 




2 3 4 

m(2(7tV )) (GeV/c 2 ) 



1 us 2 2.5 
m(n + nr\) (GeV/c ) 



FIG. 23: (a) Background-subtracted invariant mass distribu- 
tion for the four pions recoiling against the r\ in 2(7r + 7r~)?7 
events, (b) Invariant mass distribution for all tv + tt~ti com- 
bination (four entries per event) in selected 2(-k + tt~)t) can- 
didates (open histogram) and the estimated non-ISR back- 
ground (hatched). 



against the 77, after subtraction of the ISR and non-ISR 
backgrounds. The concentration around 1.5 GeV/c 2 is 
consistent with final state 77/?(1450) being one of the dom- 
inant channels in the 2(7r + 7r _ )?7 process. Figure 23(b) 
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shows the mass distribution for all neutral 7r + 7r~?7 combi- 
nations (four entries per event). Signals from the 7/(958) 
and a peak at 1.3 GeV/c 2 are evident. There are two 
candidates decaying to 777r + 7r~ and allowed by quantum 
numbers for the later: /i(1285) and 77(1295) [5]. For 
events with an entry in one of these peaks, the mass of 
the remaining 7r+7r~ pair is concentrated in the p(770) 
region, indicating that these events are predominantly 
from the 7y'(958)p and 77(1295)^ (A(1285)p). The process 
e+e~ — > /i(1285)p(770) seems to be prefered, because 
/i(1285) has the decay to 7p(770), but 77(1295) decays to 
777T7r with pions in S-wave [5] (and not well studied). We 
now study these events in detail. 



space, we obtain: 

cr = 0.18 ±0.07 nb, 
m x = 1.99 ±0.08 GeV/c 2 , 
T x = 0.31 ±0.14 GeV . 

This might be the p(2150), the next radial excitation of 
the p family, reported previously and listed in the de- 
tailed PDG tables [5]. The structure around the J/ip 
region can not be explained by the J/ip — > 7/(958)7r + 7r~ 
decay, but could be a background from other J/t/j decay 
modes with a missing ir° or undetected radiative pho- 
ton(s). 




m(it + it T|) (GeV/c 2 ) 



FIG. 24: (a) Distribution of the ir + -K~r] mass closest to the 
77' mass in the data (open histogram), along with the esti- 
mated non-ISR background contribution (hatched). The ver- 
tical lines indicate the 77' (958) signal and side band regions, 
(b) The e + e~— > 77'(958)7T + 7r _ cross section obtained via ISR. 
The line is the result of the Breit-Wigner fit described in the 
text. 



D. The 77'(958)7r+7r" Intermediate State 

To extract the contribution from the 7/(958)7r + 7r~ 
channel, we select the n + ir~ri combination with mass 
closest to the 77' (958) mass. Figure 24(a) shows the dis- 
tribution of this mass below 1.1 GeV/c 2 . A clean 77' (958) 
signal is visible. Also shown is the estimated contribu- 
tion from non-ISR background, which is small but also 
shows an 77' peak. 

We obtain a cross section in a manner similar to that 
described in Sec. IV C. We first subtract the non-ISR 
background in each mass bin, then subtract the events in 
two sidebands, 930-945 and 975-990 MeV/c 2 , from those 
in an 77' (958) signal region, 945-975 MeV/c 2 , obtaining 
a total of 120 ± 14 77'(958)"7r + 7r~ events. Repeating this 
procedure in bins of the 2(7r + 7r~)77 invariant mass and di- 
viding by the efficiency, ISR luminosity, 77—^77 branching 
fraction, and the t/'(958) — * 777r + 7r~ branching fraction of 



0.445 [5], we obtain the e 4 



■77'(958)7T+7r- 



cross sec- 



tion shown in Fig. 24(b) and listed in Table VIII. It 
shows a resonance-like behavior at around 2.1 GeV and 
a sharp drop at 3.3 GeV. Fitting a single Breit-Wigner 
function, Eq. 3 with mo = 1.5 GeV to describe the phase 



E. The /1 (1285)77+77" Intermediate State 

Figure 25(a) shows an expanded view of the 7r + 7r~?7 
invariant mass distribution [Fig. 23(b)] in the region 
around 1.3 GeV/c 2 . We fit this distribution with a 
Breit-Wigner signal function plus a second order polyno- 
mial to describe the combinatorial background, obtain- 
ing 649 ± 85 events in the peak. The fitted mass and 
width, 1.281 ±0.002 ±0.001 GeV/c 2 and 0.035 ± 0.006 ± 
0.004 GeV, are compatible with the /i(1285) parame- 
ters [5] and are not in agreement with those listed for 
the 77(1295). We conclude that contribution from the 
?7(1295)7r + 7r~ is small and more data is needed for de- 
tailed study. A similar fit to the 77' peak gives a mass 
shifted from the PDG value by -0.9 ± 0.4 MeV, from 
which we estimate a 1 MeV systematic error on the 
/1 (1285) mass. The systematic error on the width is es- 
timated by varying the shape of the combinatorial back- 
ground. 




T0.5 - 



1.1 1.2 1.3 1.4 1.5 
m(7tVr|) (GeV/c 2 ) 




E m (GeV) 



(a) Distribution of the 7r 7r 77 invariant mass in 
GeV/c 2 region for the data (points). The line rep- 



FIG. 25: 
the 1-1.5 

resents the result of the fit described in the text, with the 
shaded region representing the background component, (b) 
The e + e~ — * /i(1285)7r + 7r~ cross section obtained via ISR. 
The line is the result of the Breit-Wigner fit described in the 
text. 



We extract the number of /i(1285)7r + 7r~ events in 
80 MeV/c 2 bins of the 2(7r + 7r _ )?7 mass using similar fits 
with the /i(1285) mass and width fixed to the above val- 
ues. Due to the uncertainty in the background shape 
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TABLE VIII: Measurements of the e + e — + r/(958)7r + 7r cross section (errors are statistical only). 



(GeV) 



(nb) 



E c 



(GeV) 



(nb) 



E c 



(GeV) 



(nb) 



Ec 



(GeV) 



(nb) 



1.58 
1.66 
1.74 
1.82 
1.90 
1.98 



0.02 
0.06 
0.01 
0.07 
0.11 
0.16 



0.05 
0.03 
0.06 
0.06 
0.07 
0.06 



2.06 
2.14 
2.22 
2.30 
2.38 
2.46 



0.24 
0.10 
0.11 
-0.05 
0.03 
0.09 



0.08 
0.06 
0.05 
0.05 
0.04 
0.05 



2.54 
2.62 
2.70 
2.78 
2.86 
2.94 



0.05 
0.03 
0.03 
-0.01 
0.06 
0.04 



0.04 
0.03 
0.02 
0.02 
0.02 
0.02 



3.02 
3.10 
3.18 
3.26 
3.34 
3.42 



0.07 ± 0.03 
0.07 ± 0.03 
0.04 ± 0.02 
0.07 ± 0.02 
0.00 ± 0.02 
0.01 ± 0.01 



TABLE IX: Measurements of the e + e — * fi (1285)7r + 7r cross section (errors are statistical only). 



E c 



(GeV) 



(nb) 



E c .m. (GeV) 



(nb) 



(GeV) 



(nb) 



Ec 



(GeV) 



(nb) 



1.66 
1.74 
1.82 
1.90 
1.98 
2.06 



0.00 ± 0.00 
0.01 ± 0.11 
0.02 ± 0.18 
0.35 ± 0.20 
0.61 ± 0.23 
0.59 ± 0.24 



2.14 
2.22 
2.30 
2.38 
2.46 
2.54 



0.99 
0.89 
0.64 
0.54 
0.61 
0.25 



± 0.24 
± 0.24 
± 0.21 
± 0.17 
± 0.16 
± 0.12 



2.62 
2.70 
2.78 
2.86 
2.94 
3.02 



0.32 ± 
0.13 ± 
0.31 ± 
0.25 ± 
0.07 ± 
0.24 ± 



0.12 
0.10 
0.10 
0.08 
0.06 
0.08 



3.10 
3.18 
3.26 
3.34 
3.42 
3.50 



0.11 ± 0.08 
0.08 ± 0.05 
0.03 ± 0.04 
0.02 ± 0.03 
0.09 ± 0.05 
0.05 ± 0.04 



we assign an additional 10% systematic error on the 
number of signal events. Dividing the fitted number of 
events by the efficiency, ISR luminosity and 77—5-77 and 
/i(1285) — > rjir + ir~ branching fractions yields the cross 
section shown in Fig. 25(b) and listed in Table IX. There 
is no evidence of the J/ip decay into this mode, and the 
cross section again demonstrates resonance-like behavior 
at around 2.1 GeV. Fitting with a single Breit-Wigner 
function, Eq. 3 with mo = 1.8 GeV, we obtain: 

cr = 1.00 ±0.18 ±0.15 nb, 
m x = 2.15 ±0.04 ±0.05 GeV/c 2 , 
F, = 0.35 ±0.04 ±0.05 GeV. 

The mass and width are consistent with those obtained 
above for the 77'(958)7r + 7r _ channel, and with those listed 
in the PDG tables [5] for the p(2150), but the cross sec- 
tion is substantially larger. 

VI. THE K+K--K+-K--K FINAL STATE 

A. Final Selection and Backgrounds 

To suppress ISR 2(7r + 7r~)7r° background, we fit each 
event under that hypothesis and require xl^ ^ ^0- The 
x\k2-ktt° distribution for the remaining events is shown 
as points in Fig. 26 and the distribution for simulated 
A" + A'~7r + 7r~7r° events (open histogram) is normalized 
to the data in the region X^.k2tt-k° < 20. The hatched 
histogram represents the non-ISR background, which 
is dominated by K + K~ir + ir~ir ir a events and checked 
against data as in Sec. IV A. Since the statistics are low 
and the scale factor is consistent, we use the same value 
as for the 2(7r + 7r~)7r° final state. The largest remaining 




FIG. 26: Distribution of x 2 from the 5C fit for K + K~ -k+ -k~ 
candidates in the data (points). The open histogram is the 
distribution for simulated signal events, normalized as de- 
scribed in the text. The hatched and dashed histograms are 
the backgrounds from non-ISR events and ISR K + K ~-k + -k~ 
events, respectively, estimated as described in the text 



background is from ISR K + K~ir + ir~ events. This con- 
tribution, estimated from the simulation using the mea- 
sured cross section [11], is shown as the dashed histogram 
in Fig. 26. All other backgrounds are either negligible or 
distributed uniformly in X2K2tttt ■ We define a signal re- 
gion, X2K27ttt° <40, containing 5565 events and a control 
region, 40 < X2K27r7r° < containing 1758 events. 

Figure 27 shows the K + K~Tr + ir~7r° invariant mass 
distribution from threshold up to 4.5 GeV/c 2 for events 
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2 3 4 

m(K + K"7i + 7l"7i°) (GeV/c 2 ) 

FIG. 27: Invariant mass distribution for selected 
K + K~tv + ty~tt° candidates in the data (points). The 
hatched and open histograms represent, cumulatively, the 
non-ISR background and the ISR background from the 
control region of Fig. 26. 
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a) Invariant mass distribution for simulated 
events in the signal (open) and control 



FIG. 28 

K+K~ 7T ' 7T 7T 

(hatched) regions of Fig. 26. (b) Net reconstruction and se- 
lection efficiency as a function of mass obtained from this 
simulation; the curve represents a third order polynomial fit. 



in the signal region. The hatched histogram represents 
the non-ISR background, and the open histogram rep- 
resents the sum of all backgrounds, where the ISR-type 
background is estimated from the control region. The to- 
tal background is about 15% at low mass, but accounts 
for a large fraction of the selected events above about 
3.5 GeV/c 2 . We subtract the sum of backgrounds from 
the number of selected events in each mass bin to obtain a 
number of signal events. Considering uncertainties in the 
cross sections for the background processes, the normal- 
ization of events in the control region and the simulation 
statistics, we estimate a systematic uncertainty on the 
signal yield of less than 5% in the 1.6-3.0 GeV/c 2 region, 
but increasing to 10% in the region above 3 GeV/c 2 . 



B. Selection Efficiency 

The detection efficiency is determined in the same 
manner as in Sec. IV B. Figure 28(a) shows the simu- 
lated K + K~7r + ir~ir° invariant mass distributions in the 
signal and control regions from the phase space simu- 
lation. We divide the number of reconstructed events 
in each mass interval by the number generated in that 
interval to obtain the efficiency shown as the points in 
Fig. 28(b); the curve represents a third order polynomial 
fit to the points, which we use to calculate the cross sec- 
tion. Simulations assuming dominance of the </>7r + 7r~7r° 
and rjK + K~ channels give consistent results, and we ap- 
ply a 5% systematic uncertainty for possible model de- 
pendence, as in Sec. IV B. 

We correct for track- and 7r°-finding efficiencies, and 



the shape of the X2K2Tm° distribution as in Sec. IV B. We 
measure the kaon identification efficiency using e + e~— ► 
0(1020)7 — > K + K~-f events, as described in Ref. [11], 
and apply a correction of +(2.0 ± 2.0)% to the efficiency. 
The total efficiency correction is +8% and we estimate 
a systematic error of 10% for masses below 3.0 GeV/c 2 , 
increasing to 30% at 4.5 GeV/c 2 . 




cm. 

FIG. 29: The e + e~ — > K + K~ tv + tt~ tv° cross section as a func- 
tion of e + e~ cm. energy measured with ISR data. The errors 
are statistical only. 
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C. Cross Section for e+e 



K+K~ 



calculate the cross section for the e + e~ — > 
7r + 7r _ 7T process by dividing the number of events 
"7r + 7r _ 7r° mass bin by the corrected ef- 



We 
K+K 

in each K+K 
ficiency and differential luminosity. We show the first 
measurement of this cross section in Fig. 29, with statis- 
tical errors only, and list the results in Table X. Again, 
the energy resolution is much smaller than the bin width 
and we apply no correction. The cross section rises to 
a peak value near 1 nb at 2.5 GeV, followed by a slow 
decrease with increasing energy. The only statistically 
significant structure is the J/i/j peak. 
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FIG. 30: (a) The 7r + 7r - 7r° mass versus two-kaon-three-pion 
mass, and (b) the tv+tt~tt° mass projection for selected 
K + K~-k+tt~-k° candidates. The hatched histogram repre- 
sents the estimated non-ISR background. 
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FIG. 32: (a) The ty+tt~tv° mass distribution in the region 
below 700 MeV/c 2 for events with a K+K~ mass within 
15 MeV/c 2 of the mass (points). The histogram is the dis- 
tribution for the ISR (fyq simulation, (b) The e + e~— >0?7 cross 
section measured here in the K+ K~tv+tt~tt° final state (large 
dots), compared with the previous BABAR measurement in the 
ISR K+K -77 final state [27]. 



shows a scatter plot of the tt+tt~tt° mass versus the 
K+K~ mass in the event. A vertical band correspond- 
ing to the 0(1020) is visible, and almost all T] are pro- 
duced through the 4>rj channel, whereas the w(782) band 
is spread out across the plot. The K+K~ mass projec- 
tion in Fig. 31(b) shows a 0(1020) signal, but the non-ISR 
background distribution has no structure. 
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FIG. 31: (a) The ti + tt~ti mass versus the K + K~ mass and 
(b) the K + K~ mass projection for selected K + K~ ir + n~ it 
candidates. The hatched histogram represents the estimated 
non-ISR background. 



D. Substructure in the K+K 7r + 7r 7t° Final State 



E. The 4>r\ Intermediate State 

Requiring the K + K~ mass to be within ±15 MeV/c 2 
of the nominal 0(1020) mass and plotting the mass of the 
recoiling 7r + 7r~7r° system below 700 MeV/c 2 , we obtain 
the distribution shown in Fig. 32(a). The distribution 
from the ISR 4>rj simulation is also shown, and this chan- 
nel can account for all entries below 700 MeV/c 2 . Count- 
ing events with a three-pion mass in the 0.5-0.6 GeV/c 2 
region in bins of the K + K~ir + ir~ir° mass, and divid- 
ing by the corrected efficiency [Fig. 28(b)], differential 
luminosity and — > K + K~ and rj — > ir+ir~ir° branching 
fractions, we obtain the e + e~— >0?7 cross section shown 
in Fig. 32(b) and listed in Table XI. 

The cross section shows a rise from threshold to a peak 
value of about 2 nb at around 1.7 GeV, followed by a 
monotonic decrease with increasing energy. This mea- 
surement is consistent with the more precise BABAR mea- 
surement in the K + K~~jj final state [27], which is also 
shown in Fig. 32(b) 



Figure 30(a) shows a scatter plot of the ir + ir~ir 
mass versus the K + K~ir + 7T~ir° mass in the selected 
K + K~7T + ir~7T° events. There are horizontal bands cor- 
responding to the 77 signal and w(782) resonance. The 
7r + 7T~7r° mass projection, Fig. 30(b), shows ij and u> 
peaks, as well as a small signal for the 0(1020). The non- 
ISR background contribution is shown as the hatched his- 
togram and also contains 7/ and lj signals. Figure 31(a) 



F. The u>(782)K+K~ Intermediate State 

Figure 33(a) shows the ir + 7T~ir° mass distribution 
in the region near the lj(782) mass for all selected 
K + K~7T + ir~ir° events in the data and the estimated 
non-ISR contribution. We first subtract the non-ISR 
background in each mass bin, then subtract the events 
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TABLE X: Measurements of the e + e — > K + K tt + -k it cross section (errors are statistical only). 



£c.m. (GeV) 


a 


(nb) 


E c . m . (GeV) 


a 


(nb) 


E e . m . (GeV) 


a 


(nb) 


Bc.m. (GeV) 


a 


(nb) 


1.6125 


0.02 


± 


0.04 


2.3375 


0.89 


± 


0.13 


3.0625 


0.86 


± 


0.15 


3.7875 


0.38 


± 0.10 


1.6375 


0.08 


± 


0.06 


2.3625 


0.73 


± 


0.13 


3.0875 


7.01 


± 


0.34 


3.8125 


0.18 


± 0.08 


1.6625 


0.13 


± 


0.08 


2.3875 


0.82 


± 


0.13 


3.1125 


4.31 


± 


0.28 


3.8375 


0.36 


± 0.09 


1.6875 


0.27 


± 


0.10 


2.4125 


0.75 


± 


0.13 


3.1375 


0.83 


± 


0.15 


3.8625 


0.36 


± 0.08 


1.7125 


0.17 


± 


0.08 


2.4375 


0.85 


± 


0.13 


3.1625 


0.92 


± 


0.15 


3.8875 


0.20 


± 0.08 


1.7375 


0.09 


± 


0.06 


2.4625 


1.04 


± 


0.13 


3.1875 


0.75 


± 


0.13 


3.9125 


0.37 


± 0.08 


1.7625 


0.32 


± 


0.10 


2.4875 


1.02 


± 


0.14 


3.2125 


0.57 


± 


0.13 


3.9375 


0.20 


± 0.07 


1.7875 


0.16 


± 


0.07 


2.5125 


0.95 


± 


0.14 


3.2375 


0.62 


± 


0.14 


3.9625 


0.33 


± 0.08 


1.8125 


0.34 


± 


0.09 


2.5375 


0.85 


± 


0.15 


3.2625 


0.47 


± 


0.12 


3.9875 


0.18 


± 0.08 


1.8375 


0.36 


± 


0.09 


2.5625 


0.82 


± 


0.14 


3.2875 


0.43 


± 


0.13 


4.0125 


0.11 


± 0.07 


1.8625 


0.32 


± 


0.09 


2.5875 


0.63 


± 


0.12 


3.3125 


0.50 


± 


0.13 


4.0375 


0.14 


± 0.07 


1.8875 


0.20 


± 


0.09 


2.6125 


1.20 


± 


0.15 


3.3375 


0.55 


± 


0.13 


4.0625 


0.39 


± 0.08 


1.9125 


0.54 


± 


0.11 


2.6375 


0.73 


± 


0.13 


3.3625 


0.42 


± 


0.12 


4.0875 


0.13 


± 0.07 


1.9375 


0.54 


± 


0.11 


2.6625 


0.71 


± 


0.12 


3.3875 


0.48 


± 


0.11 


4.1125 


0.16 


± 0.07 


1.9625 


0.63 


± 


0.11 


2.6875 


0.93 


± 


0.14 


3.4125 


0.37 


± 


0.11 


4.1375 


0.19 


± 0.07 


1.9875 


0.57 


± 


0.11 


2.7125 


0.73 


± 


0.13 


3.4375 


0.50 


± 


0.12 


4.1625 


0.13 


± 0.07 


2.0125 


0.39 


± 


0.10 


2.7375 


0.80 


± 


0.13 


3.4625 


0.29 


± 


0.11 


4.1875 


0.13 


± 0.08 


2.0375 


0.49 


± 


0.11 


2.7625 


0.77 


± 


0.13 


3.4875 


0.52 


± 


0.12 


4.2125 


0.14 


± 0.06 


2.0625 


0.71 


± 


0.13 


2.7875 


0.64 


± 


0.14 


3.5125 


0.45 


± 


0.11 


4.2375 


0.21 


± 0.06 


2.0875 


0.56 


± 


0.11 


2.8125 


0.71 


± 


0.13 


3.5375 


0.42 


± 


0.10 


4.2625 


0.04 


± 0.08 


2.1125 


0.59 


± 


0.12 


2.8375 


0.82 


± 


0.13 


3.5625 


0.42 


± 


0.10 


4.2875 


0.10 


± 0.06 


2.1375 


0.59 


± 


0.12 


2.8625 


0.64 


± 


0.13 


3.5875 


0.53 


± 


0.11 


4.3125 


0.20 


± 0.07 


2.1625 


0.66 


± 


0.12 


2.8875 


0.83 


± 


0.14 


3.6125 


0.18 


± 


0.09 


4.3375 


0.14 


± 0.07 


2.1875 


0.60 


± 


0.13 


2.9125 


0.65 


± 


0.13 


3.6375 


0.25 


± 


0.11 


4.3625 


0.28 


± 0.07 


2.2125 


0.48 


± 


0.11 


2.9375 


0.76 


± 


0.14 


3.6625 


0.34 


± 


0.10 


4.3875 


0.05 


± 0.06 


2.2375 


0.61 


± 


0.12 


2.9625 


0.63 


± 


0.13 


3.6875 


0.60 


± 


0.12 


4.4125 


0.25 


± 0.08 


2.2625 


0.49 


± 


0.11 


2.9875 


0.93 


± 


0.15 


3.7125 


0.43 


± 


0.10 


4.4375 


0.41 


± 0.09 


2.2875 


0.84 


± 


0.13 


3.0125 


0.78 


± 


0.14 


3.7375 


0.36 


± 


0.10 


4.4625 


0.19 


± 0.06 


2.3125 


0.55 


± 


0.12 


3.0375 


0.62 


± 


0.13 


3.7625 


0.38 


± 


0.10 


4.4875 


0.10 


± 0.07 



TABLE XI: Measurements of the e + e — > <^(1020)?7 cross section (errors are statistical only). 



(GeV) 



(nb) 



E c 



(GeV) 



(nb) 



E c 



(GeV) 



(nb) 



E c 



(GeV) 



(nb) 



1.56 
1.60 
1.64 
1.68 
1.72 
1.76 
1.80 



0.00 ± 0.41 



0.57 
0.70 
2.28 
1.53 
2.02 
1.23 



0.33 
0.35 
0.61 
0.48 
0.54 
0.41 



1.84 
1.88 
1.92 
1.96 
2.00 
2.04 
2.08 



0.65 
0.63 
0.85 
0.58 
0.57 
0.22 
0.54 



0.29 
0.28 
0.32 
0.26 
0.25 
0.16 
0.24 



2.12 
2.16 
2.20 
2.24 
2.28 
2.32 
2.36 



0.63 
0.31 
0.30 
0.00 
0.20 
0.10 
0.29 



0.26 
0.18 
0.18 
0.20 
0.14 
0.10 
0.17 



2.40 
2.44 
2.48 
2.52 
2.56 
2.60 
2.64 



0.10 ± 0.10 
0.19 ± 0.13 
0.00 ± 0.19 
0.28 ± 0.16 
0.00 ± 0.19 
0.19 ± 0.13 
0.09 ± 0.09 



in two sidebands, 690-735 and 825-870 MeV/c 2 , from 
those in an u> signal region, 735-825 MeV/c 2 . The mass 
distribution of the recoiling K + K~ pair after this back- 
ground subtraction shows no resonant structure. Di- 
viding by the corrected efficiency, differential luminosity 
and lo — > 7r + 7T~7r° branching fraction, we obtain the first 
measurement of the e + e~ -^uj (782) K + K~ cross section, 
shown in Fig. 33(b) and listed in Table XII. 

The cross section rises from threshold to a peak value 
of about 0.55 nb at about 2 GeV, then decreases with 
increasing energy except for peaks at the J/ip and ip(2S) 
masses. The events in the latter peak are partly due to 
the i>(2S)-> J/tPn+n', J/iP-tK+K-n decay. 



VII. THE K + K-n + n-ri FINAL STATE 

A. Final selection and backgrounds 

To suppress ISR 2(7r + 7r~)?7 background, we fit each 
event under that hypothesis and require xl™ > 30. The 
X2K2TT7] distribution for the remaining events is shown as 
points in Fig. 34, and the distribution for simulated ISR 
4>t]^ K + K~ir + ir~ir° events (open histogram) is normal- 
ized to the data in the region X2K2mr° < ^0- We do not 
simulate K + K~n + n~ri events, but we expect the resolu- 
tion and efficiency to be indistinguishable from the gen- 
erated mode. The hatched histogram represents the non- 
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TABLE XII: Measurements of the e+e" -> cj(782)if + if" cross section (errors are statistical only). 



E c . m . (GeV) 


a 


(nb) 


E c . m . (GeV) 


a 


(nb) 


E c .m. (GeV) 


a 


(nb) 


E c .m. (GeV) 


a 


(nb) 


1.80 


0.09 


± 


0.04 


2.48 


0.31 


± 


0.07 


3.16 


0.07 


± 


0.04 


3.84 


0.04 


± 0.02 


1.84 


0.33 


± 


0.08 


2.52 


0.27 


± 


0.06 


3.20 


0.12 


± 


0.04 


3.88 


0.01 


± 0.01 


1.88 


0.25 


± 


0.06 


2.56 


0.18 


± 


0.05 


3.24 


0.17 


± 


0.05 


3.92 


0.01 


± 0.01 


1.92 


0.52 


± 


0.09 


2.60 


0.20 


± 


0.05 


3.28 


0.06 


± 


0.04 


3.96 


0.01 


± 0.01 


1.96 


0.60 


± 


0.10 


2.64 


0.20 


± 


0.05 


3.32 


0.07 


± 


0.04 


4.00 


0.02 


± 0.01 


2.00 


0.45 


± 


0.09 


2.68 


0.18 


± 


0.06 


3.36 


0.06 


± 


0.03 


4.04 


0.02 


± 0.01 


2.04 


0.43 


± 


0.09 


2.72 


0.10 


± 


0.04 


3.40 


0.01 


± 


0.02 


4.08 


0.00 


± 0.01 


2.08 


0.61 


± 


0.10 


2.76 


0.16 


± 


0.05 


3.44 


0.07 


± 


0.03 


4.12 


0.00 


± 0.01 


2.12 


0.37 


± 


0.08 


2.80 


0.18 


± 


0.05 


3.48 


0.04 


± 


0.02 


4.16 


0.00 


± 0.01 


2.16 


0.28 


± 


0.07 


2.84 


0.12 


± 


0.05 


3.52 


0.07 


± 


0.03 


4.20 


0.02 


± 0.01 


2.20 


0.39 


± 


0.08 


2.88 


0.12 


± 


0.05 


3.56 


0.08 


± 


0.03 


4.24 


0.00 


± 0.01 


2.24 


0.36 


± 


0.08 


2.92 


0.17 


± 


0.05 


3.60 


0.07 


± 


0.03 


4.28 


0.00 


± 0.01 


2.28 


0.17 


± 


0.06 


2.96 


0.17 


± 


0.05 


3.64 


0.05 


± 


0.03 


4.32 


0.04 


± 0.02 


2.32 


0.22 


± 


0.06 


3.00 


0.16 


± 


0.05 


3.68 


0.14 


± 


0.04 


4.36 


0.02 


± 0.01 


2.36 


0.34 


± 


0.07 


3.04 


0.05 


± 


0.04 


3.72 


0.09 


± 


0.03 


4.40 


0.02 


± 0.01 


2.40 


0.23 


± 


0.06 


3.08 


0.51 


± 


0.09 


3.76 


0.02 


± 


0.02 


4.44 


0.00 


± 0.01 


2.44 


0.19 


± 


0.06 


3.12 


0.24 


± 


0.06 


3.80 


0.00 


± 


0.02 


4.48 


0.00 


± 0.01 




0.7 0.8 0.9 

m(7[ + 7i"jt°) (GeV/c 2 ) 



(GeV) 



FIG. 33: (a) Expanded view of the 7r + 7r~7r° invariant mass 
distribution [Fig. 30(b)] in the region near the oj(782) for all 
selected events in the data (open histogram) and the esti- 
mated non-ISR background (hatched). The vertical lines de- 
limit the u> signal region and sidebands, (b) The e + e~ — > 
u>(782)K + K ~ cross section obtained via ISR. 



ISR background, which is dominated by K + K~7r + ir~7r a ri 
events and is evaluated from the simulation using the 
same scale factor as for the K + K~ Tt + -K~ n a final state. 
We define a signal region, "X^ K2 . KV < 40, containing 375 
events and a control region, ^0<X2K2ttt 1 < ^' containing 
162 events. We subtract the non-ISR background and the 
ISR-type background, estimated from the control region, 
to obtain a number of signal events. 



B. Cross Section for K + K ' 



We calculate the cross section for the process e + e~- 
K + K~ir + ir~ri as described in Sec. IV C, by dividing 
the number of background-subtracted events in each 
K + K^TT + -K^r] mass bin by the corrected detection ef- 
ficiency and differential luminosity. Since the model de- 
pendence of the acceptance is small, we use the efficiency 




100 



X (2K2jcti) 



FIG. 34: Distribution of X 2 from the 5C fit for K + K~ ^ ^ t] 
candidates in the data (points). The open histogram is the 
distribution for simulated signal events, and the hatched his- 
togram is the estimated background from non-ISR events. 



for K + K~ir + ir~ir° events shown in Fig. 28(b), divided 
by the r\ — > 77 branching fraction and with an increased 
systematic error of 5%. 

We show the cross section as a function of energy in 
Fig. 35 with statistical errors only. This is the first mea- 
surement of this cross section, which shows a rise from 
threshold to a maximum value of about 0.2 nb at about 
2.8 GeV, followed by a monotonic decrease with increas- 
ing energy, except for a prominent peak at the J/ip mass. 
The systematic errors are similar to those for the other 
modes presented here, totalling about 10% at all energies. 
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FIG. 35: The e + e~— * K + K~ n + n~ r) cross section as a func- 
tion of cm. energy measured with ISR data. Only statistical 
errors are shown. 
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FIG. 36: (a) The 7r + 7r _ 77 mass versus the K + K~ tt + ty~ r] mass, 
and (b) the ir + -K~r] mass projection for selected K + K~ ix + ix~ -q 
candidates. The hatched histogram represents the estimated 
non-ISR background. 
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FIG. 37: (a) The 7r + 7r _ 7j mass versus the K + K~ mass and 
(b) the K + K~ mass projection for selected K + K~ -k + it~ r) 
candidates. The hatched histogram represents the estimated 
non-ISR background. 
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FIG. 38: Raw invariant mass distribution for all selected 
e + e~ — > 2(7T + 7r _ )7r events in the charmonium region. The 
line represents the result of the fit described in the text. 



C. Substructure in the K+K-TT+TT-rj Final State 

Figure 36(a) shows a scatter plot of the 7r + 7T~?7 mass 
versus the K + K~-K + -K~ri mass and Fig. 36(b) shows the 
7T + 7r~?7 mass projection. A horizontal band and peak, re- 
spectively, corresponding to the 77' (958) are visible. The 
non-ISR background, shown as the shaded histogram in 
Fig. 36(b), is small, but may include a few 77'. Fig- 
ure 37(a) shows a scatter plot of the tt + tt^t] mass versus 
the K + K~ mass in the event. A vertical band corre- 
sponding to the 0(1020) is visible, and almost all 77' (958) 
are produced through the 4>rf channel. The K + K~ mass 
projection in Fig. 37(b) shows a 0(1020) signal, but the 
non-ISR background distribution has no resonant struc- 
ture. Due to the low statistics, we do not study the mass 
dependence of these channels. 



VIII. THE CHARMONIUM REGION 

The data at masses above 3 GeV/c 2 can be used to 
measure or set limits on the branching fractions of nar- 
row resonances, such as charmonia, and the narrow J/ip 
and i/}(2S) peaks allow measurements of our mass scale 
and resolution. Figures 38, 39, 40 and 41 show the in- 
variant mass distributions for the selected 2(ir + ir~)ir Q , 
2(7T + 7r~)?7, K + K~ir + ir~ir° and K + K~ir + ir~i] events, 
respectively, in this region, with finer binning than in 
the corresponding Figs. 2, 21, 27 and 35. We do not 
subtract any background, since it is small and nearly uni- 
formly distributed. Signals from the J/ip and i/)(2S) are 
visible in all four distributions. 

We fit these distributions using a sum of two Gaussian 
functions to describe each of the J/ip and ip(2S) signals 
plus a polynomial to describe the remainder of the dis- 
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TABLE XIII: Measurements of the e + e" -> if 7r + 7r rj cross section (errors are statistical only). 



£cm. (GeV) 


a 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


E c .m. (GeV) 


a 


(nb) 


£c.m. (GeV) 


a 


(nb) 


2.1125 


0.00 


± 


0.00 


2.7125 


0.09 


± 


0.06 


3.3125 


0.11 


± 


0.08 


3.9125 


0.08 


± 0.08 


2.1375 


0.08 


± 


0.06 


2.7375 


0.12 


± 


0.07 


3.3375 


-0.08 


± 


0.09 


3.9375 


-0.03 


± 0.06 


2.1625 


0.08 


± 


0.05 


2.7625 


0.14 


± 


0.08 


3.3625 


0.18 


± 


0.09 


3.9625 


0.04 


± 0.06 


2.1875 


0.00 


± 


0.00 


2.7875 


0.23 


± 


0.11 


3.3875 


0.12 


± 


0.09 


3.9875 


-0.01 


± 0.06 


2.2125 


0.04 


± 


0.04 


2.8125 


0.17 


± 


0.09 


3.4125 


0.00 


± 


0.05 


4.0125 


0.08 


± 0.06 


2.2375 


0.04 


± 


0.04 


2.8375 


0.06 


± 


0.08 


3.4375 


0.00 


± 


0.10 


4.0375 


0.03 


± 0.07 


2.2625 


-0.02 


± 


0.02 


2.8625 


0.05 


± 


0.07 


3.4625 


0.18 


± 


0.08 


4.0625 


0.11 


± 0.07 


2.2875 


0.04 


± 


0.04 


2.8875 


0.07 


± 


0.07 


3.4875 


0.17 


± 


0.10 


4.0875 


0.05 


± 0.07 


2.3125 


0.00 


± 


0.00 


2.9125 


0.03 


± 


0.08 


3.5125 


0.05 


± 


0.05 


4.1125 


0.05 


± 0.06 


2.3375 


0.00 


± 


0.00 


2.9375 


0.22 


± 


0.12 


3.5375 


0.13 


± 


0.11 


4.1375 


-0.03 


± 0.05 


2.3625 


-0.02 


± 


0.02 


2.9625 


0.06 


± 


0.10 


3.5625 


0.01 


± 


0.08 


4.1625 


-0.07 


± 0.05 


2.3875 


0.07 


± 


0.05 


2.9875 


0.34 


± 


0.13 


3.5875 


-0.01 


± 


0.08 


4.1875 


0.00 


± 0.06 


2.4125 


-0.02 


± 


0.05 


3.0125 


0.10 


± 


0.09 


3.6125 


0.09 


± 


0.09 


4.2125 


0.00 


± 0.05 


2.4375 


0.00 


± 


0.00 


3.0375 


0.16 


± 


0.10 


3.6375 


0.09 


± 


0.10 


4.2375 


0.10 


± 0.08 


2.4625 


0.17 


± 


0.09 


3.0625 


0.25 


± 


0.13 


3.6625 


0.10 


± 


0.07 


4.2625 


0.06 


± 0.04 


2.4875 


0.02 


± 


0.04 


3.0875 


1.63 


± 


0.26 


3.6875 


0.18 


± 


0.10 


4.2875 


0.12 


± 0.08 


2.5125 


0.09 


± 


0.06 


3.1125 


0.95 


± 


0.22 


3.7125 


0.21 


± 


0.09 


4.3125 


0.04 


± 0.06 


2.5375 


0.08 


± 


0.08 


3.1375 


0.24 


± 


0.10 


3.7375 


0.06 


± 


0.05 


4.3375 


0.15 


± 0.07 


2.5625 


0.09 


± 


0.06 


3.1625 


0.20 


± 


0.11 


3.7625 


-0.01 


± 


0.08 


4.3625 


0.02 


± 0.07 


2.5875 


0.16 


± 


0.08 


3.1875 


0.05 


± 


0.08 


3.7875 


0.00 


± 


0.07 


4.3875 


0.08 


± 0.07 


2.6125 


0.05 


± 


0.05 


3.2125 


0.36 


± 


0.14 


3.8125 


0.08 


± 


0.08 


4.4125 


0.06 


± 0.07 


2.6375 


0.26 


± 


0.10 


3.2375 


0.17 


± 


0.13 


3.8375 


-0.04 


± 


0.09 


4.4375 


-0.02 


± 0.08 


2.6625 


0.14 


± 


0.07 


3.2625 


0.20 


± 


0.11 


3.8625 


0.15 


± 


0.07 


4.4625 


0.03 


± 0.06 


2.6875 


-0.04 


± 


0.06 


3.2875 


0.04 


± 


0.09 


3.8875 


0.15 


± 


0.07 


4.4875 


0.04 


± 0.11 
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m(2(7tV)Ti) (GeV/c 2 ) 

FIG. 39: Raw invariant mass distribution for all selected 
e + e~ — ► 2(Tv + Tv~)ri events in the charmonium region. The 
line represents the result of the fit described in the text. 




m(K + K"7tV/) (GeV/c 2 ) 

FIG. 40: Raw invariant mass distribution for all selected 
e + e~ — + K + K~ tt + it~ tt° events in the charmonium region. 
The line represents the result of the fit described in the text. 



tribution. We fix the two Gaussian combination shape 
by fitting simulated events, but let the overall mean and 
width float in the fit, along with the amplitude and the 
coefficients of the polynomial. In all cases, the fitted 
mean values are within 1 MeV/c 2 of the PDG [5] J/tp 
and tp(2S) masses, and the widths are consistent with 



the simulated resolutions within 10%. 

In the J /if) peak the fits yield 4990±79 2(tt+tt-)t: 
events, 85±14 2(tt + tt~)ti events, 768±31 K+K^n+TT-n 
events and 72.9±9.4 K + K~n + -K~ri events. From the 
number of observed events in each final state /, Nj/^f, 
we calculate the product of the J /ip branching fraction 
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Dividing by the world average value of T 



V>(2S) 



[5], 



3.25 



3.5 



3.75 



m(K + K"7i;Vr|) (GeV/c ) 



FIG. 41: Raw invariant mass distribution for all selected 
e + e~ —> K + K~ it + -k~ events in the charmonium region. The 
line represents the result of the ht described in the text. 



to / and the J/ip electronic width: 



&J/ip-+f ' ^ie — 



N 



6tt 2 ■ dC/dE ■ tf(mjj^) ■ C 



(5) 



where dC/dE = 65.6 ± 2.0 nb _1 /MeV and e/(m J/v ,) 
are the ISR luminosity and corrected selection efficiency, 
respectively, at the J/tp mass and C is a conversion 
constant. We estimate efficiencies of 0.105 for the 
2(7r + 7T~)7r° and 2(ir + ir~)r] final states, and 0.046 for the 
K + K~ir + ir~ir° and K + K~ir + ir~i] modes, with system- 
atic uncertainties of about 5% for modes with ir° and 
about 7% for modes with -q. Adding the error on the ISR 
luminosity in quadrature, we assign a 6% (7.5% for rf) 
overall systematic uncertainty on each product. 

Using re/ V ' = 5.55± 0.14 keV [5], we obtain the branch- 
ing fractions listed in Table XIV, along with the mea- 
sured products and the current PDG values. The sys- 
tematic errors include the 2.5% uncertainty on vil^ . 
The J/tp — ► K + K~Tr + ir~rj channel has not been pre- 
viously observed and the K + K~ir + ir~ir° and 2(ir + ir~)ri 
branching fractions are consistent and competitive with 
the PDG values. However, we find a 2(7r + 7r _ )7r° branch- 
ing fraction 4.8 standard deviations higher than the PDG 
value. 

In the ip(2S) peak the fits yield 410±30 
2(ir + TT~)ir° events, 15.6±7.6 2(7r + 7r - )?7 events, 
31.8±11.9 K+K-tt+it-tt events and 7.0±4.0 
K + K~ir + ir~ri events. Using a calculation analo- 
gous to Eq. 5, with dC/dE = 84.0 ± 2.5 nb~7MeV 
and e(ip(2S)) = 0.0965 (0.0400) for the 2(7r+7r-)vr 
and 2(7r + 7r~)?7 (K + K~Tr + ir~ir° and K + K~n + ir~ri) 
modes, we obtain the product of the ip(2S) branching 
fractions to these final states and its electronic width. 



we 

obtain the branching fractions listed in Table XIV. The 
2(7r + 7r~)?7 and K + K~Tt + Tr~r) branching fractions are 
first measurements, and the K + K ~Tr + ir~ir° branching 
fraction is consistent with the PDG value. 

However, we find a 2(7r + 7r~)7r° branching fraction 7.6 
standard deviations higher than the PDG value. We 
note that some of the observed ip(2S) could be due to 
the decay chain ip(2S) — > J/ipir + Tr~, J/ip — ► ir + ir~ir°, 
and we use this chain to check our result. The scat- 
ter plot of the 7r + 7r~7r° mass closest to the J/ip mass 
versus the 2(7T + 7r _ )7r° mass in Fig. 42(a) shows a clus- 
ter corresponding to this decay chain. We select events 
with a three-pion mass within 50 MeV/c 2 of the J/tp 
mass (lines in Fig. 42(a)) and plot their 2(tt + tt^)tt° mass 
in Fig. 42(b). A fit yields 256±17 ip(2S) events, and 
using the well measured ip(2S) — > J/ipir + ir~ branch- 
ing fraction of 0.318±0.06 [5], we calculate a J/ip — > 
7T + 7r~7r° branching fraction Bj/^^^+^-^o — (2.36 ± 
0.16 ± 0.16)% that is consistent with our previous mea- 
surement -Bj/^-tTr+ir-TT = (2.19 ± 0.19)% [8] as well as 
with the current PDG value. We obtain significantly 
higher values for both the J/ip and ip(2S) branching 
fractions to 2(7r + 7r~)7r° compared to previous experi- 
ments [5]. A similar difference was reported for the 
J/ip — * 7T + 7r~7r° decay in recent experiments [8, 28]. 
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FIG. 42: (a) The tv + ty tv° mass closest to J/ip mass versus the 
five-pion mass for selected 2(n + n' -)n° events, (b) The five- 
pion mass for events with a three-pion mass within 50 MeV 
of the J/ip mass. 
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FIG. 43: Background-subtracted 2(7r + 7r - )7r° invariant mass 
distributions in the charmonium region for events with a 
7T + 7r _ 7r° mass in the (a) r\ and (b) w(782) mass region. 
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TABLE XIV: Measurements of the J/tp and tp(2S) branching fractions. 
Measured Measured J/tp or ip(2S) Branching Fraction (10~ 3 ) 





Quantity 


Value (eV) 


Calculated, this work 


PDG2006 


r J/4> 

J- ee 


^J/i/>^2(7r+ir-)ir° 


303 


±5 ±18 


54 


6 ±0.9 ±3.4 


33.7 ±2.6 


r J/i> 

J- ee 




47 


8 ±3.1 ±3.2 


9. 


7 ±0.6 ±0.6 


7.2 ±1.0 


J- ee 


B.J/il>~^r)Tr+TV- ' Bn^3ir 





51± 0.22± 0.03 





40± 0.17± 0.03 


0. 193± 0.023 


r J/i> 

1 ee 




5 


16± 0.85± 0.39 


2 


35± 0.39± 0.20 


2.26 ±0.28 


r J/i> 

1 ee 


Bj/ii^K+K-TT+n-TT 


107 


±4.3 ±6.4 


19 


2 ±0.8 ±1.5 


12.0 ±3.0 


r J/« 

J- ee 


Bj/ip->4>ri ■ B<t,^K+K- ' £>I7^3tt 





84± 0.37± 0.05 


1. 


4 ±0.6 ±0.1 


0.74 ±0.08 


1 ee 


^J/ip—tuiK+K- ' £>w->3tt 


3 


3 ± 1.3 ±0.2 


1. 


36±0.50±0.10 


1.9 ±0.4 


r J/i> 

J- ee 


Bj/lp-*K+ K-TT+TT-J) ' By^^y 


10 


2 ± 1.3 ±0.8 


4. 


7 ±0.6 ±0.4 


no entry 


P*(2S) 
-L ee 


^i/>(2S)-*2(7r + 7r-)7r° 


29 


7 ± 2.2 ± 1.8 


12 


±0.9 ±0.7 


2.66 ±0.29 


r *(2S) 
-L ee 


^i/'(2S)^J/!/ J 7r+7r- ' t3j/i>-^:iTT 


18 


6 ± 1.2 ± 1.1 


23 


6 ± 1.6 ± 1.6 


20. 2 ± 1.4 


pV(2S) 




2 


69±0.73±0.16 


1. 


22± 0.33± 0.07 


0.66 ±0.17 


p*(2S) 


^V(2S)-»J/i/'i7 ' ^7->37r ■ Sj/^-^+j,- 


1 


11±0.33±0.07 


33 


4 ±9.9 ±2.0 


30.9 ±0.8 


P*(2S) 
-L ee 


^V(2S)^2(7r+7i--)?) ' Bri^ll 


1 


13± 0.55± 0.08 


1. 


2 ±0.6 ±0.1 


no entry 


p*(2S) 


"i> (2S)-»Jf + if - 7T+ 7T - 7T° 


4 


4 ± 1.3 ±0.3 


1. 


8 ±0.5 ±0.1 


1.24 ±0.10 


r *(2S) 
1 ee 


^V(2S)-*if+K-7r+7T-i7 ' Br/^iy 


1 


2 ±0.7 ±0.1 


1. 


3 ±0.7 ±0.1 


no entry 



We are also able to measure J/tp and tp(2S) branching 
fractions for some of the submodes studied above. Fig- 
ure 43 shows expanded views of the 2(ir + ir~)n mass 
distribution in the charmonium region for the rjir+Tr^ 
and ujir + ir~ intermediate states. Our fits yield 8.9±3.8 
and 788±51 J/tp decays to rjiT + ir~ and ujir + n~ , respec- 
tively, and 14.2±4.2 and 37±10 tp(2S) decays. We list 
the corresponding products and branching fractions in 
Table XIV. The tp(2S) — ► rj-K+Tr^ branching fraction is 
very small [5] and the observed events are from the decay 
chain tp(2S) — > J/tpr], J/tp — > , rj — > 7r + 7r~7r°. The 

result in Table XIV assumes this decay chain. We also 
observe 6.0±2.7 and 24±9 events in the J/tp peaks for the 
<jyq and K + K ~u modes, respectively. The corresponding 
products are also listed in Table XIV. 



IX. SUMMARY 

The photon energy and charged particle momentum 
resolutions together with the particle identification capa- 
bilities of the BABAR detector permit the reconstruction 
of e + e"->- 2(7T + 7r-)7r°, 2(ir + iT-)r], K+K-Tr+n-n and 
K + K ~7r + ir~ri events produced at low effective e + e~ cm. 
energy via ISR in data taken in the Y(4S) mass region. 
Luminosity and efficiency can be understood with 6-10% 
accuracy, so that ISR production yields useful measure- 
ments of R, the ratio of the hadronic to di-muon cross 
section values, used for the (g— 2) p calculations. 

Our measurements of the e + e~ — > 2{tt + tt~)tt cross 
section represent a significant improvement upon exist- 
ing data in both energy range and precision. In addi- 
tion, these data provide new information on hadron spec- 
troscopy. The observed e + e~ '— > lutt^ ~tt~ and t]'k + 'k~ cross 
sections show evidence of resonant structures around 1.4- 



1.7 GeV/c 2 , which were previously observed by DM2 and 
interpreted as o>(1420) and w(1650) resonances. We ob- 
tain new measurements of the parameters of these reso- 
nances, which confirm the results of our previous study 
of ISR 7T + 7r~7r events, that the o>(1650) is substantially 
narrower than currently listed in PDG. 

We also use this final state to make the first measure- 
ments of the e + e~^ u/o(980) and p(770)37r cross sec- 
tions. In the latter events, there is an iso- vector resonant 
structure with m = 1.243 ± 0.012 ± 0.020 GcV/c 2 and 
T = 0.410 ± 0.031 ± 0.030 GcV in the thrce-pion system 
recoiling against the p(770). 

We present the first measurements of the cross sections 
for e + e-->2(7r+7r-)r?, rj' (958)t: + n~ and /i(1285)7r + 7r _ . 
We measure the mass and width of the /i(1285), and 
observe a candidate p(2150) resonance with m p ( 2 i 5 o) = 
2.15 ± 0.04 ± 0.05 GeV/c 2 and r p(2150 ) = 0.35 ± 0.04 ± 
0.05 GeV. 

We present the first measurements of the cross sections 
for e + e~ — > K + K~ 7r + ir~ tt° , K + K~u> and K + K~Tr + Tr~r/. 
Using the latter final state, we measure a contribution 
from e + e~^(prj consistent with our measurement in the 
K+K-- n final state [27]. 

The final states analyzed in this paper, based on 
232 flT 1 of BABAR data in the 1.0-4.5 GeV/c 2 mass 
range, are already better in quality and precision than 
the direct measurements from the DCI and ADONE ma- 
chines, and do not suffer from the relative normalization 
uncertainties which seem to exist for direct measurements 
of these final states. 

The ISR events also allow a study of J/tp and tp(2S) 
production, and the measurement of fifteen products of 
branching fractions into observed modes and the e+e~ 
width of the J/tp or tp(2S). Three of these are first mea- 
surements, and two others are the most precise measure- 
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merits to date, thanks to our relatively small systematic 
error due to acceptance. We observe substantial dis- 
crepancies with respect to the previous experiments in 
the J/tp -> 2(n + TT-)TT° and ip(2S) -> 2(7r+7r-)7r° decay 
modes. 
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